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1.  INTRODUCTION 


The  US  Army  has  had  a  continuing  interest  in  the  design  of  spinning 
projectiles  with  liquid  payloads.  Many  of  these  developmental  shell  have 
shown  dramatic  instabilities  in  their  pitching  and  yawing  motion.  Initially, 
these  instabilities  have  been  identified  by  large  range  losses  incurred  dur¬ 
ing  firing  trials.  In  1962  Karpov  made  direct  angular  motion  measurements 
of  liquid-payload-induced  instabilities  in  a  20mm  projectile  fired  in  BRL's 
Aerodynamics  Range1.  In  1973  Mark  and  Mermagen2  used  instrumented  solar 
sensors  and  telemetry  units  to  observe  liquid-payload-induced  instabilities 
in  155mm  shell.  Since  that  time  all  developmental  shell  with  liquid  payloads 
have  teen  tested  with  sunsonde  instrumentation  and  a  variety  of  strange  be¬ 
haviors  has  been  observed.  A  complete  listing  of  ail  yawsonde  data  reports 
is  given  in  the  Bibliography  at  the  end  of  this  report. 

In  1959  Stewartson  published  a  theoret  cal  paper  on  the  stability  of  a 
spinning  liquid-filled  top3.  This  paper  assumed  a  right  circular  container 
partially  or  fully  filled  with  an  inviscid  fluid.  The  liquid  was  assumed  to 
be  fully  spun  up  and  in  steady  state  motion.  This  motion  was  assumed  to  be 
a  circular  or  spiral  motion  at  a  frequency  set  by  the  top's  static  moment 
and  spin  rate.  The  Stewartson  theory  predicted  liquid  eigenfrequencies  that 
were  to  be  avoided  in  order  to  have  stable  angular  motion  of  the  top.  Ac¬ 
cording  to  the  theory,  liquid  moments  would  become  infinite  for  coning  motion 
at  any  of  the  eigenfrequencies. 

In  1965  Karpov1*  made  additional  20mm  firings.  All  shell  in  this  series 
had  the  same  fast  frequencies  but  the  payload  eigenfrequencies  were  varied 
by  use  of  different  cavity  fineness  ratios.  A  resonance  undamping  rate  was 
observed  but  was  of  a  much  smaller  amplitude  and  at  a  slightly  lesser  fre¬ 
quency  than  that  predicted  by  Stewartson. 

At  about  that  time  Wedemeyer5  introduced  a  boundary  layer  modification 
to  the  Stewartson  theory  and  computed  complex  liquid  eigenvalues.  Since  the 
Aerodynamics  Range  flights  were  too  short  to  allow  the  liquid  to  be  fully 
spun  up  and  in  steady  state  coning  motion,  Karpov  developed  the  use  of  a 
free  liquid-filled  gyroscope  to  measure  yaw  growth  rates  near  resonance.  He 
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found  an  exceptionally  good  agreement  with  Wedemeyer  values  at  a  Reynolds 
number  of  520,000  and  fair  agreement  at  a  Reynolds  number  of  520Q6. 

The  success  of  Karpov's  gyroscope  experiments  led  to  an  extensive  use 
of  this  technique.  A  complete  listing  of  gyroscope  data  reports  is  also 
included  in  the  Bibliography.  This  excellent  experimental  work  had  the  un¬ 
fortunate  effect  of  biasing  most  of  the  later  theoretical  and  experimental 
work  toward  understanding  liquid-filled  gyroscopes  and  the  application 
to  projectiles  was  treated  as  a  side  effect. 

A  second  difficulty  with  the  later  gyroscope-oriented  work  was  a  tunnel 
vision  concentration  on  liquid  eigenfrequencies.  This  was  caused  by  the 
great  success  of  Wedemeyer' s  modification  of  Stewartson's  inviscid  eigen¬ 
frequencies.  The  basic  aim  of  any  liquid  payload  theory  should  be  the  calcu¬ 
lation  of  the  complete  moment  the  liqu’d  payload  exerts  on  the  pitching  and 
yawing  projectile  in  flight.  Wedemeyer's  complex  eigenfrequencies  identify 
frequency  and  damping  rate  pairs  for  which  the  liquid  pressure  is  infinite. 
For  coning  motion  near  any  of  these  pairs  the  liquid  moment  is  primarily  due 
to  the  pressure  at  the  edge  of  the  boundary  layer  and  is  dominated  by  a  sim¬ 
ple  pole  function.  This  pole  is  an  excellent  approximation  at  high  Reynolds 
number,  but  at  lower  Reynolds  number  it  becomes  quite  poor  even  though  the 
boundary  layer  assumptions  are  still  valid.  The  pressure  at  the  edge  of  the 
boundary  layer  has  to  be  computed  without  the  pole  approximation.  In  addi¬ 
tion  an  increment  in  pressure  through  the  rotating  boundary  layer  on  the 
lateral  wall  must  be  computed,  as  well  as  the  shear  on  both  the  lateral  and 
end  walls. 

It  is  the  aim,  then,  of  this  report  to  give  the  general  formulation  of 
the  effect  of  liquid  payload  motion  on  projectile  stability  and  to  compute 
the  liquid  moments,  pressures  and  wall  shears  for  small-ampl itude  liquid  mo¬ 
tion  with  boundary  layers  but  without  the  unnecessary  mathematical  approxima¬ 
tions  of  the  Stewartson-Wedemeyer  theory.  The  results  of  this  improved 
Stewartson-Wedemeyer  theory  will  be  compared  with  all  available  published 
gyroscope  data  for  Reynolds  numbers  down  to  as  low  as  2400.  Moreover,  the 
theory  will  be  extended  to  the  special  case  of  a  fully-filled  cylinder  with 
a  central  rod7'*.  Finally,  a  survey  of  extensions  of  the  theory  to  partially 
spun-up  liquids  and  other  special  cases  will  be  given. 
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2.  PROJECTILE  DYNAMICS 


Two  coordinate  systems,  both  of  which 
axis  of  symmetry,  are  commonly  used:  the 


have  X-axes  along  the  projectile  s 
missile-fixed  XYZ  system  and  the 


aeroball istic  XYZ  non-rolling  system  with  the  Z-axis  initially  downward.  IT 

we  introduce  earth-fixed  axes  X  Y  Z  with  the  X  -axis  initially  along  the 

e  e  e  e 

velocity  vector  and  Zg  downward,  a  unit  vector  along  the  positive  X-axis  has 
earth-fixed  components  (n^,  n^,  n^).  The  angle  of  attack  a  in  the  non¬ 


rotating  system  is  the  angle  in  the  XZ  plane  from  the  X-axis  to  the  velocity 
vector,  and  the  angle  of  sideslip  8  is  the  angle  in  the  XY  plane  from  the 
X-axis  to  the  velocity  vector.  Thus  for  a  straight  trajectory  and  small 
angles,  these  angles  are  the  negatives  of  the  direction  cosines  nZE  and  nYp 
respectively. 


The  primary  lateral  force  on  the  projectile  is  the  normal  force,  whic/ 
can  be  easily  expressed  in  terms  of  complex  variables10: 

Fy  +  i  =  -  (1/2)  pV*S  CN  l 

where 

K  =  6  +  i  a 

and  where  the  other  symbols  are  defined  in  the  List  of  Symbols. 

For  an  approximately  straight  trajectory,  the  usual  linear  aerodynamic 
moment  can  be  expressed  as  the  sum  of  three  terms: 


(2.1) 

(2-2) 


♦  i  »  U/2)  ,>VJSi  jju./V)  CHp  -  iCH  ]  l 

a 


-  i  (CM  ♦  CMJ  U  i/V) 

q  ‘x 

The  first  tern  is  the  very  important  Magnus  moment,  which  is  a  viscous  side 
moment  caused  by  the  spin  and  the  angle  of  attack.  The  second  tem  is  the 
static  moment  which  for  most  projectiles  causes  an  increased  angle  in  the 
plane  of  the  total  angle  of  attack.  The  last  tem  is  the  damping  rodent, 
which  usually  resists  the  angular  velocity.  For  simplicity  we  will  Reelect 
the  effect  of  drag  on  the  angular  notion  and  assume  V  to  be  constant  in 
Equations  (2.1)  and  (2.3). 


.Vujtc*, 


if." 


■  rat 

,  a: 


Hr 


:  ’sme  ‘ 


Ft  >»  /• 

.  ;•  «*  (  3  v  *• 
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For  small  angles  the  usual  dynamics10  yield  the  following  differential 
equation  for  £  in  terms  of  an  arbitrary  force  and  moment: 


y  • 


(I  d/dt  -  iU)(F;  +  IF;) 

- sv^ — - 


(2.4) 


For  the  force  and  moment  of  Eqs.  (2.1)  and  (2.3)  this  reduces  to 


£  +  (H  -  io<j»)£  -  (M  +  io<(>T)£  =  0 


(2.5) 


where 


H  *  (pSi/2m) 


CN  '  ky  (CM  +  W 
a  J  q  a 


(v/O 


(PS£3/2I  )  CM  (V/£)2 

j  a 

r-  ^  o 


M 

T  =  (PS£/2m) 


CN  +  k*  CM 
a  pa 


(V/£) 


h/ly 


The  solution  to  Equation  (2.5)  is  an  epicycle  which  generates  the  angular 
motion  as  the  sum  of  two  rotating  and  damping  or  undamping  two-dimensional 
vectors : 


5 


i+2 

+  K2e 


(2.6) 


where 


in  (Kj/ltj  )  -  .j  tj  |{|  t 


b.  =  A.  +  r.*t 

vjo  y 


j  -  (a/2)  1  ±  ^l-(l/sg)j  Ssg  =  d2  ^2/4M 


t  jH  -  ctT 

**  (2tj  -  uhjiH 


Note  that  for  coning  motion  in  the  direction  of  spin,  Tj->0  while  e.>0  for  in¬ 
creasing  K.  and  e..<0  for  decreasing  K..  For  coning  motion  in  the  direction 
opposite  to  the  spin,  the  inequalities  are  reversed. 
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In  analyzing  the  effect  of  a  moving  internal  component  on  the  angular 
motion  of  a  spinning  projectile11,  we  found  it  convenient  to  consider  only 
that  part  of  the  moment  exerted  by  the  internal  component  at  one  cf  the  two 
frequencies  of  the  projectile's  angular  motion.  Tor  .teadv-state  linear 
liquid  motion,  this  part  will  be  the  total  liquid  moment.  For  nonsteady  or 
nonlinear  liquid  motion,  this  part  will  consist  of  two  average  components  of 
the  actual  liquid  moment.  If  we  now  non-dimensional ize  this  liquid  moment 

by  the  liquid  mass  m^ ,  the  spin  rate  <*•  and  the  maximum  liquid  container  di¬ 
ameter  2a,  the  following  expression  for  the  liquid  moment  can  be  obtained: 


MLY  +  iMLZ  =  mLa2^2 


lc 


t  C,M  K  e 

l  LM  i 
l 


1  +  \CLM 


K  e 
2 


(2.7) 


For  linear  fluid  motion,  C,  M  should  depend  on  x  ,  e  ,  time,  Reynolds 

LMi  1  i 

number,  fill  ratio,  the  shape  of  the  cavity,  and  the  direction  of  the  spin.  A 
similar  remark  applies  to  C^M  .  The  Xj  's  appear  explicitly  in  definition 

(2.7)  since  the  moment  should  vanish  for  Xj  =  0. 


It  should  be  noted  that  the  C(  „  are  complex  quantities  whose  imaginary 

L  J 

parts  represent  in-plane  moments  causing  rotation  in  the  plane  of  exp(i<j>.) 

J 

and  whose  real  parts  represent  side  moments  causing  rotations  out  of  the 
plane  of  exr.  «<• We,  therefore,  introduce  the  following  definition  for  the 

real  and  imaginary  parts  of  C^M  and  explicitly  express  the  effect  of  the 

0 

direction  of  spin10: 


CLM  •  =  yCLSM.  +  1  CL IM . 
J  J  J 


(2.8) 


where  and  are  real  and  represent  the  liquid  side  moment  and 

1  J  .  * 

liquid  in-plane  moment  contributions,  respectively,  and  where  c 

The  special  values  of  these  coefficients  for  infinitely  viscous  or 
frozen  liquid  can  be  obtained  from  Equation  (2.4)  with  the  external  moments 
neglected. 


<Iy  +  ILy>C  -  lid,  ♦  IlxH  -  0  (2.9) 

11.  C.  H.  tfurphy ,  "Influence  of  Mooing  Internal  Farts  on  Angulo}'  Motion  of 
Spinning  Projectiles,"  Journal  of  Gunianco  and  Control,  Vol.  1,  Marck- 
April  1978,  pp.  117-122.  (See  also  PEI,  Memoi'andwn  Rcpoi't  2731  dated 
February  1977.  AD  037338.) 
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'  <npppu  i- " 


where  1^  and  I^x  are  transverse  and  axial  moments  of  inertia  for  the  frozen 

liquid.  Comparing  Equations  (2.4)  and  (2.9)  we  see  that  if  the  terms  involv¬ 
ing  the  frozen  liquid  moments  of  inertia  are  taken  to  the  right  side,  they 
can  be  identified  as  i  times  the  liquid  moment. 


HLY  +  iMLZ  -  -« V  • 


(2.10) 


The  epicylic  solution  of  Equation  (2.6)  with  | e j j «1  can  now  be  used  to  pro¬ 
vide  frozen  liquid  values  of  the  liquid  side  moment  and  liquid  in-plane 
moment  coefficients. 


CLIM. 

J 

CLSM . 

J 


.  ‘  2Tj!tyl 

itijT7 


(2.11 ) 

(2.12) 


For  a  circular  cylinder  of  length  2c  and  center  located  a  distance  h  forward 
of  the  projectile's  center  of  mass 


I  .  c2  +  3ii 2 

^ +  ___5 — 


(ej/2)|l 


(2.13) 

(2.14) 


A  simple  interpretation  for  e  follows  from  the  observation  that  for 
moderate  damping  2kg  is  approximately  the  fractional  change  in  i(.  in  one 

>J 

cycle.  If  we  restrict  this  change  to  be  less  than  20",  e  should  be  less 
than  .03,  and  the  frozen  would  be  less  than  .015. 

In  general,  however,  the  liquid  moment  of  Eq.  (2.7)  should  be  combined 
with  the  aerodynamic  force  and  moment  of  Eqs.  (2.1)  and  (2.3)  to  give  a  some¬ 
what  more  complicated  differential  equation  for  £: 


C  +  (H  -  io,$)£  -  (M  +  io$T)e 


(2.15) 
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If  the  epicycle  solution  of  Eq.  (2.6)  is  substituted  in  Eq.  (2.15),  new  re¬ 
lations  for  frequency  and  damping  can  be  obtained. 


tj  =  (a/2) 


fj  -  H)J  Vfj!  -  <Vsg) 


where 


fj  ■  1  +  (V'Vm 


arid 


a  ** 


£j  =  ' 


x,  (H  +  H,  )  -  cT 

J _ h _ 

(2xj  -  a) tj  III 


where 


Hl  =  -<pS£3/2Iy )  CLMq  (V/A) 
J  J 


'LMq  .  =  °i  H.SM . 
3 

2m^a2 |lj 


L  pSi‘ 


(2.16) 


(2.17) 


As  can  be  seen  from  Eq.  (2.17),  the  liquid  side  moment  has  the  same  effect  on 
the  damping  of  the  angulai  motion  as  the  aerodynamic  damping  moment.  The 
coefficient  is  introduced  so  cha'  the  relative  size  of  the  aerodynamic 


damping  moment  and  the  liquid  side  moment  can  be  directly  evaluated.  The 
direct  impact  of  the  liquid  side  moment  on  the  dampinq  per  cycle  can  be  seen 

from  Eq.  (2.17)  for  H  =  T  =  3. 


cj  s  (mj_a2/Ix) (2x^0  -  1  )“l  Cj 


LSM 


j 


(-.18) 


For  the  fast  mode  the  coefficient  of  CLSM  is  positive  and  a  positive  side 
moment  causer  an  undamping  of  this  motion  Similarly  a  negative  C^SM  will 

undamp  the  slow  mode.  As  we  shall  see,  the  linear  lie*  id  motion  theory 
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usually  yields  a  positive  side  moment,  and  thus  only  the  fast  mode  motion 
is  adversely  affected  by  the  liquid  side  moment. 


Table  1.  Parameter  Values  for  Five  Army  Projectiles 


Payload 

Projectile 

.. 

Diameter 

(mm) 

c/a 

a 

mLaZ 

I 

X 

°L 

White 

M416 

105 

2.67 

.17 

.36 

350 

Phosphorus 

M328 

107 

2.82 

.11 

.41 

230 

Re  =  4-40  x  106 

XM825 

155 

4.60 

.08 

.12 

150 

Binary 

Chemical 

M687 

155 

4.52 

.08 

.07 

80 

Re  =  1-7  x  106 

XM736 

203 

. 

3.98 

.12 

.09 

90 

In  Table  1  the  parameters  c/a,  o,  m^a2/^  and  are  given  for  five 

Army  projectiles.  The  first  three  are  smoke  projectiles  containing  white 
phosphorus,  which  is  liquid  for  temperatures  above  110°F,  and  the  remaining 
two  have  special  liquid  payloads.  Since  t  lies  between  c/2  and  o,  we  see 

that  the  range  of  interest  is  .04  to  .17.  m^a2/Ix  and  oL  are  much 
larger  for  the  smoke  shell  due  to  HP's  greater  density. 

Equation  (2.18)  can  be  used  to  determine  a  lower  bound  on  the  side 
moment  coefficient  corresponding  to  a  significant  yaw  growth  rate  of  13£ 
per  cycle,  i.e.,  e  =  .02.  For  t/c  =  3/4  and  m.a2/I  =  .08,  this  lower  bound 
on  is  .125.  For  m^a2/Ix  =  .40,  which  is  appropriate  to  the  older  WP 

shell,  this  lower  bound  is  .025.  Thus  our  theoretical  prediction  of  this 
liquid  side  monrent  coefficient  should  at  least  attempt  to  achieve  an  accura¬ 
cy  better  than  half  the  lower  value,  i.e.,  errors  less  than  .01. 


3.  EQUATIONS  OF  LIQUID  MOTION 

We  will  consider  a  projectile  with  a  cylindrical  payload  cavity  with 
radius,  a,  and  height,  2c.  The  axis  of  the  cylinder  is  collinear  with  the 
projectile's  axis,  and  its  center  is  located  a  distance,  h,  from  the  pro¬ 
jectile's  center  of  mass.  If  the  cavity  is  partially  filled,  the  liquid  is 
fully  spun  up,  and  the  centrifugal  force  is  large  compared  to  the  aerodynamic 
forces,  the  liquid  will  fill  the  space  between  the  outer  cylindrical  wall  ard 


16 


an  inner  cylindrical  free  surface  with  radius,  b.  The  ratio  of  the  volume 
of  this  inner  cylinder  to  the  volume  of  the  complete  payload  cavity  is 
b2/a2.  The  fill  ratio  for  the  payload  cavity  is,  therefore,  l-b2/a2  and  will 
be  denoted  by  f.  However,  m,  will  always  be  the  liquid  mass  in  a  fully- 
filled  cavity. 

The  objective  of  the  linear  theory  is  to  predict  the  liquid  moment  re¬ 
sponse  to  coning  or  spiral  motion  of  the  form 


Ha 

K  =  K.e  J  j  =  1  or  2  (3.1) 

J 


=  Kes* 


where 


5  =  +  0  Tj 

« 

(j)  =  (j)t 


The  vector  between  the  center  of  mass  of  the  projectile  md  any  other 
points  on  the  projectile  can  be  described  in  the  aeroballlstic  cylindrical 

coordinates  of  this  vector  by  (x,  r,  e).  Cartesian  coordinates  of  this 
vector  in  the  earth-fixed  coordinates  would  then  be  (x  .  y  .  z  1.  For 

e  c 

simplicity  we  will  omit  the  subscript  V*  for  these  earth-fixed  coordinates 
throughout  this  report.  Relations  between  the  earth-fixed  Cartesian  coordi¬ 
nates  and  the  aeroballistic  cylindrical  coordinates  take  on  quite  simple 
forms  for  small  K,. 

J 

x  =  x  +  rK .  cos  ($.  -  8)  (3.2) 

J  J 


y  =  r  cos  8  -  K.x  cos  $.  (3.3) 

s? 


z  =  r  sin  6  -  K.x  sin  (3.4) 

V  V# 
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If  cylindrical  coordinates  with  respect  to  the  earth-fixed  axes  are 
denoted  by  (x,  r,  e),  the  following  simple  relations  between  the  two  sets  of 
polar  coordinates  follow  from  Eqs.  (3.3  -  3.4)  for  small  K ^ . 

**  «K» 

r  =  r  -  KjX  cos  -  e)  (3.5) 

sin  (8-8)=  (xK./r)  sin  (a.  -  e)  (3.6) 

J  J 

The  cylindrical  components  of  the  velocity  of  any  point  on  the  projectile  in 

t  *  •  • 

aeroballistic  coordinates  are  x  =  r  =  0,  re  =  r<t>.  In  earth-fixed  coordinates 
they  can  be  obtained  by  differentiating*  Eqs.  (3.2,  3.5  -  3.6). 

Vx  =  R|l(s  -  i)rK  es*  ’  19  j  (3.7) 

Vr  =  -r{*(s  -  i)xK  es^  "  18  }  (3.8) 

VQ  -  k  +  R  {U(s  -  i)xK  es*  '  18  }  (3.9) 


where 


rI  \  = 


n 


is  the  real  part  of  a  complex  quantity. 


We  will  nov#  make  the  very  restrictive  assumption  that  the  liquid  is  in 
steady-state  response  to  the  coning  and  spinning  motion  of  the  projectile. 
Theoretical  studies  12 *  13  have  been  made  and  are  in  progress  to  determine 
the  effect  of  partially  spun-up  liquid,  and  an  experimental  study  of  the 
transient  response  to  coning  motion  has  been  made14.  These  studies  show 
that  spin-up  and  cone-up  effects  are  large  and  important  to  a  complete  under¬ 
standing  of  the  liquid  payload  stability  problem. 

Nevertheless,  we  will  assume  that  the  liquid  velocity  components  and 
liquid  pressure  have  the  same  dependency  on  time  and  0  as  the  velocity  com¬ 
ponents  of  points  on  the  projectile  and  introduce  four  small  dimensionless 
functions  of  r  and  x:  u  ,  v  ,  w c  and  p  . 

due  b 


'See  Appendix  A  Jor  details. 


12.  M.  Lynn,  "Free  Oscillations  of  a  Liquid  Puring  Spin-Up t"  BRL  Report 
1665,  August  1$?5 .  AD  769716, 


13. 


C.  U.  Kitchens,  Ur.,  h\ 
in  a  Rotating  Cylinder: 
Report  ARBRL - TR- 0 2 6 8 1 , 


Gerber,  and  R.  Sedney,  "Cecil lai 
Part  I.  Solid-Body  Rotation," 
June  1978.  AD  AGS77S8. 


tens  of  a  ; 
BRL  Technic 


14.  U.  F.  D’Amico,  W.  G.  Seims,  and  T.  !L  Regers,  "Pressure  Measurements  of 
a  Rotating  Liquid  for  Impulsive  Coning  Motion "  BRL  Menotxtndum  Report 
in  publication.  (See  also  /SJAA  Paper  82-6249,  Jan  1982. ) 
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=  4 


v„  =  R  (u_es^  *  r  J  (a*) 


Vr  =  R{VS^1e}<4) 


V0  =  r$  +  R  |wses*  "  ie  ]>  (a$) 


P  = 


PL 


£  +  R{p§es^ 


16  }  (pj_a2<^s 


(3.10) 

(3.11) 

(3.12) 

(3.13) 


Eqs.  (3.10  -  3.13)  can  now  be  placed  in  th-  linearized  unsteady  Navier- 
Stokes  equations  and  the  continuity  equation  yield 


3p  .  , 

(s  -  1)v$  -  2ws  +  a  -J*  =  yRe 


azve,  2a2iwc 

v9%-V  +  -^S 


r‘  J 


(3.14) 


iap 

(s-i)ws.2vs-  — 


=  yRe 


ci 


a2wr  2azi v_ 

*,  2W _ I _ §. 

v8  w  2  ra 


(3.15) 


3Ps 

3X 


J 

(s  -  i)u  +  a—  =  yRe'1  vn2u 


e  “s 


3(rv  )  3u 

~Jr~  -  iws  +  r  FT  “  0 


(3.16) 

(3.17) 


where 


v  2  =  a2 
o 


U 

[_  ar 


■  3  . 

irT  rar  3x2 


■a 


4.  80UNDARY  LAYER  SOLUTION 


Wederoeyer5  made  the  assumption  that  the  velocity  components  and  pressure 
could  each  be  expressed  as  the  sum  of  inviscid  and  viscous  terms.  The  invis- 
cid  terms  satisfy  Eqs.  (3.14  -  3.17)  for  Re"1  -  0  over  the  entire  cylinder 

except  for  a  small  ooundary  layer  region  near  the  cylinder  walls,  while  the 
viscous  terns  satisfy  the  boundary  layer  versions  of  Eqs.  (3,14  -  3.17). 
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Although  Wedemeyer  considers  the  effect  of  these  boundary  layer  terms  only 
on  the  liquid  eigenvalues,  this  report  will  consider  all  their  contributions 
to  the  liauid  moment.  Since  the  effect  of  negative  spin  can  easily  be  found 
from  Eq.  (2.8),  we  will  only  consider  positive  spin  (y  =  1)  for  the  remainder 
of  this  report. 


=  usi 

+ 

usv 

(4.1) 

=  V  • 

+ 

V„„ 

(4.2) 

SI 

sv 

=  w  . 

+ 

w 

(4.3) 

SI 

sv 

=  Psi 

+ 

psv 

(4.4) 

On  the  lateral  wall  r  =  a,  then,  the  usual  boundary  layer  approximations* 
reduce  Eqs.  (3.14  -  3.17)  to: 


•  -3*  ■  2wsv 

(4.5) 

(s  -  Owsv  -  a^Re 

(4.6) 

(s  -  i)u  =  a»e_1 
v  *  sv  ar2 

(4.7) 

3(rv  ,)  3u 

sv  *  sv 

- — -  =  iw  -  r 

ar  sv  ax 

(4.8) 

Far  from  the  lateral  wall,  u$,  vs>  ws,  p$  must  vanish.  At  the  wall  the 

velocities  must  be  those  required  by  Eqs.  (3.7  -  3.9).  The  viscous  tan¬ 
gential  velocities  can  be  determined  and  a  condition  for  the  inviscid  normal 
velocity  obtained.* 

wsv  =  [^d  +  is)(x/a)K  -  ws.je  (r'a}/a6a  (4.9) 

usv  =  _[(i  -  s)l<  +  (4.10) 

*  See  Appendix  8  for  details. 
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For  r  =  a, 


v  .  -  as 


8vsi 


si  a  3r 


(i  -  s){x/a)K 


(Ml) 


where 


a 


^2(1  +  is) 


At  the  end  walls,  x  =  x  ~  -  =  ±1 , 


3^V 

(s  -  i)vcw  -  2wcw  =  a2Re_1  - 

'  SV  SV  3X2 


(4.12) 


-i  32wsv 


(s  -  *>wsv  +  2vsv  =  a2Re  -IP 


(4.13) 


3psv 

a  — —  =  0 
3X 


(4.14) 


3usv  . 
r  ~W  1Wsv 


3(rVsv) 

3T~ 


(4.15) 


Once  again  the  solution  for  the  tangential  viscous  velocities  and  a  relation 
for  the  normal  inviscid  velocity  can  be  obtained. 


\v  *  ivsv  a  -  (wsi  +  ivsi)e 


-a(l  +  x) 


(4.16) 


W  -  1  V 
SV  SV 


•  -  [wsi  -  ivsj  -  2(1  +  isltM-2)  K  ]  e'9°  *  X)  (4J7) 


For  x  =  ±1 , 


3usi 


tU  +  5  — r-  3  -  (i  -  s)(r/a)K 
3x 


(4.18) 
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where 


a  =  (c/a)  <sa*  virTTiiTrrMiT 

B  =  i (c/a)  6a'  yftf  -  i s )/( 1  +  is) 


<$  =  _(a/c)6a 
2-^1  +  is 


1  -  is 
^  3  +  is 


i 


It  is  interesting  to  note  that  according  to  Eq.  (4.14)  the  usual 
boundary  condition  of  no  pressure  change  through  the  boundary  layer  is  pres¬ 
ent.  Eq.  (4.5)  shows  that  this  is  not  the  case  on  the  lateral  wall.  The 
pressure  at  the  wall  differs  from  that  at  the  edge  of  the  boundary  layer  by 
Psv(a).  This  pressure  difference  can  be  computed  by  inserting  w sv  as  given 

by  Eq.  (4.9)  in  Eq.  (4.5)  and  integrating. 


Psv(a,x)  =  2«a  [  (1  *  i s) (c/a) (x  +  h/c)  K  -  w$i  (a,x) 


(4.19) 


5.  INVISCID  SOLUTION 

The  inviscid  terms  are  solutions  of  Equations  (3.14  -  3.17)  for  Re"1  =  a 

These  four  equations  can  be  easily  manipulated  to  yield  a  partial  differential 
equation  for  p^  and  three  equations  for  the  three  velocity  components  in 

terns  of  psi- 


(s  -  i)2 


32P 


SI 


3rJ 


3P 


si 

r3r 


_s1 

r2 


(s2  -  2  is  + 


32p  .. 
3)  -  51 
1  3x2 


(5.1) 


(s  -  i)u 


st 


3p 


SI 


3X 


-a 


(5.2) 


(5.3) 


3Pc-j  2iap  • 

(s2  -  21s  +  3)v$i  =  -(s  -  0  a  -^  +  — ^ 

,  3Psi  ai(s  -  i)ps1 

(s2  -  2is  +  3)w$i  =  2a  + - - - 21  (5.4) 


On  the  outer  walls  of  the  container  the  boundary  conditions  are  given  by 
Equations  (4.18)  and  (4.11).  These  equations  can  be  rewritten  by  use  of 
Equations  (5.2)  and  (5.3). 


For  x  =  ±1 


3P«H  -  32P<H 

__|i  +  *  — s2i  *  -fe  -  1)2  (c/a)  (r/a)  K 
3x  u  3x2 


(5.5) 


For  r  =  a 


i  (1  ♦  «a)Psi  -  [s  -  i(l  -  2  8a)]  3%!-. 


2i 


a26a(s  -  i) 


9?psi 

Ir7-  =  -  (sz  +  1)(s  -  3i)(x/a)  K 


(5.6) 


On  the  inner  free  boundary,  r  =  b,  the  pressure  must  be  a  constant. 


Se  =  i|*  v  a  +  v  ^  .  v  =  o 

dt  3t  x  3x  r  3r  a  r3e 


(5.7) 


or 


(s  -  i)  Psi  +  (r/a)  v$.  «  0  (5.8) 

Note  that  for  a  fully  filled  projectile.  Equation  (5.8)  requires  p^  to  be 

zero  for  r  =  b  =  0.  Equation  (5.8)  can  now  be  simplified  by  use  of  Equation 
(5.3). 

For  r  =  b 

[(s2  +  l)(s-3i)  +  2i] psi  -  (s  -  i)  r  =  0  (5.9) 
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An  obvious  solution  to  Eq.  (5.1)  which  satisfies  Eq.  (5.5)  is 

Psi  =  -  (s  -  i)2(xr/a2)  K  (5. 10) 

We,  therefore,  assume  the  general  solution  to  be 

Psi  =  -(c/a)[(s  -  i)2(x/c)(r/a)  +  SRk(r)Xk(x)J  K  (5.11) 

Substitution  of  Eq.  (5.11)  in  Eq.  (5.1)  shows  that  Xk  is  a  linear 
combination  of  a  sine  and  a  cosine.  Eq.  (5.5)  can  be  used  to  completely 


specify  this 

combination*  for  small  6Q. 

Xk  =  cos  (lax),  k  even 

(5.12) 

=  sin  (kxx),  k  odd 

(5.13) 

=1,  k  =  0 

(5.14) 

where. 

, 

X  =  (tt/2)[1  +  6C1 

Corresponding  to  these  X^'s,  Eq.  (5.1)  gives  the  general  form  of  the  R^’s. 

Eqs.  (5.6)  and  (5.9)  can  then  be  used  to  completely  specify  the  R^'s.  In 

order  to  do  this,  x  in  Eqs.  (5.6)  and  (5.11)  must  be  replaced  by  a  series  in 
the  X^'s.  This  is  easy  when  the  X^'s  are  orthogonal.  Unfortunately,  for 

6c  not  equal  to  zero  the  X^'s  of  Eqs.  (5.12  -  5.13)  are  not  orthogonal.  We 

can,  however,  approximate  x  by  a  least  squares  fit  to  a  truncated  series  in 

V 

N 

Vk**1  (5-,5) 

k  «  0 


*See  Appendix  C  for  details*.  As  car.  be  even  there,  our  solution 

is  a  special  fc rm  cf  X  F  that  satisfies  the  inhonogenecus  for*  of  To.  ($.$'■ 

In  general,  inhomogeneous  boundary  conditions  are  rot  so  friendly,  and  much 
ihore  algebraic  labor  is  required. 
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where 


a  =  h/c 
o 

a^  =  0  for  k  even, 

and  ak  is  computed  in  Appendix  D  for  k  odd.  A  first  approximation  for 
small  6C  is  the  usual  orthogonal  relation 

\ '  Vbkk  (5-16) 

where 

bk  *J  x  Xk  (x)  dx 


bkk  “  h  \  dx 

p^.  now  assumes  the  slightly  simpler  form 

N 

psi  =  ‘(c/a)  *  K(>':  +  (S  ‘  1)2  (r/a)2k]  (5J7) 

k  s  0  J 

Equation  (5.1)  can  now  be  used  to  get  the  general  form*  of  the  R^'s. 

R.  3  (h/c)  j\  r/a  ♦  F^a/rj  (5.18) 

For  k  odd 

-  ak  j^Ekd  ( k ). r / c )  ♦  FfcY  (k\r/c)J  (5.19) 


teiz:  -s. 


25 


where 


s2  -  21  s  +  3 
.  (s-i)2 


X2 


^k’  Rk  sre  coeff  icients  to  be  determined  by  boundary  conditions^, 

J  is  a  Bessel  function  of  the  first  kind  of  order  n, 
n 

Y  is  a  Bessel  function  of  the  second  kind  of  order  n. 
n 


The  radial  functions  R^(r)  must  satisfy  boundary  conditions  (5.6)  and  (5.9). 

Direct  substitution  of  Eq  (5.17)  in  these  equations  yields  the  following 
conditions: 

2i  (1  f  3a)  Rk(a)  -  [s-i  (1  -  26#)]  a  R'k(a)  (5.20) 

+  (s  -  l)6aa2  R"k(a)  =  2aks  (s  -  i)(s  -  31) 

[(s2  +  l)(s  -  31)  +  21 ]  Rk(b)  -  (s  -  i)  b  R 1 k ( b)  (5.21) 

=  -  ak  (b/a)  s2  (s  -  i)2  (s  -  3i) 


Equations  (5.20  -  5.21)  can  be  used  to  find  E  and  F  for  specific  values  of 

0  0 

(c/a),  (b/a),  and  s.  For  non-zero  k,  R^(r)  is  a  sum  of  Bessel  functions  with 
derivatives  given  by  the  following  equations15: 

rR‘k  =  (kxr/c)ak  [  E^kxr/c)  +  FkY#  (kXr/c)]-  RR  (5.22) 

r?R"k  =  Rk  [l  "  (^r'c)2]  -  rR'k  (5.23) 

With  these  relations  the  boundary  conditions  (5.20  -  5.21)  yield  pairs  of 
linear  equations  for  Ek  and  Fk  which  can  be  quickly  solved. 


15.  N.  W.  MoLaohlan ,  Bessel  Functions  ■  for  Engineers Oxford  University  Press } 
London ,  1955. 
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5.  PRESSURE  MOMENT 

The  major  components  of  the  liquid  moment  are  due  to  the  pressure  on 
the  lateral  and  end  walls  of  the  container.  Lesser  components  are  due  to 
the  viscous  wall  shear  on  the  lateral  and  end  walls.  Thus  the  liquid  moment 
coefficient  can  be  given  as  a  sum  of  four  terms. 


tC1m  =  rn  „  +  nr  +  nv.  +  in  (6.1) 

LM  pa  pe  va  ve  \u*w 


The  sum  of  the  first  two  terms  is  the  pressure  moment  coefficient,  m  ,  and 

p 

will  be  computed  in  this  section.  The  wall  shear  moment  coefficient  will  be 
computed  in  the  next  section. 

By  use  of  Eq.  (3.5),  we  can  express  the  fluctuating  part  of  the  inviscid 
pressure  given  by  Eq.  (3.13)  as 


Jr  «j[Psi-(^)K]es*-^| 

(6.2) 

=  R  jcpeHp  K  '  ie)  J 


Equation  (5.17)  can  now  be  used  to  give  the  following  series  for  the  pressure 
coefficient: 


10  ” 

Cpe  =  -<c/a>]C  Xk  (x">  K<r)  +  <s  -  21  >s  (r/a)akl 
k  =  0  1  J 

(6.3) 

■  -(c/a)  ^  Xk  M  Cpk  (r) 


The  pressure  moment  coefficient  on  the  lateral  wall  can  be  computed  by 
an  integral  of  the  real  pressure  over  this  wall,  with  the  appropriate  lever 
arm. 
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p  £ 


X  2U 

=  i  (c/a ) (2tt  K)"1  e”s<^  J  j  (x  +  h/c)eie  r{ 

-X  0 


s$  -  ie . 

V  } 


=  i (c/2a) 


/  M 
J  -1 


dx  +  (h/c)2  m, 


r  =  a 


pjth 


do  dx 


r  *  a 


(6.4) 


where 

Vh  =  1(c2/2ah)  J  CP 


„  <a>  dx"i  V  =  Ve'*P  +  psv 


Since  the  complex  pressure  dependence  on  x  is  a  sum  of  sine  functions,  this 
integral  can  be  easily  evaluated. 


The  end  wall  pressure  moment  coefficient  is  the  difference  of  two  simi¬ 
lar  integrals  on  each  end  wall. 


m  =  -i  (a/c)(2itK)  a  3  e 

PC 


-ix  =  1 


f  a/*11Tie  /  -  s*  +  i(o 

j  e  R{cPKe 


P  Vzd9  dr 


*  -1(a/2c)  a 


-  3 


/ 


e  p  C. 


b  0 
x  =  1 

x  =  -1 


X  = 


r2  dr 


(6.5) 


This  moment  coefficient  involves  the  integrals  of  Bessel  functions  but  these 
particular  integrals  can  be  easily  obtained  in  closed  form.15 

Stewartson3  incorrectly  used  a  complex  ap  in  his  pressure 

integral  and,  therefore,  his  moment  calculation  lacks  the  (j)  factors  of 
Eqs.  (6.4  -  6.5).  Later  he  made  a  similar  error  of  a  factor  of  two  in  com¬ 
puting  a  complex  direction  cosine  so  that  these  cancelling  errors  give  the 
correct  yaw  growth  rates,  e-.  These  yaw  growth  rates  when  modified  by  Wede- 

meyer5  gave  outstanding  agreement  with  gyroscope  measurements  for  large 
Reynolds  numbers. 


In  order  to  compute  the  pressure  moment  coefficient,  it  is  necessary  to 
determine  the  parameters  Ek  and  Fk  from  the  boundary  condition  equations 

(5.20  -  5.21).  For  Stewartson's  inviscid  case  (<$a  =  0)  and  constant  ampli¬ 
tude  coning  motion  (e  =  0)  the  conditions  for  the  coefficients  of  the  Bessel 
functions  (k  f  0)  are* 


c  +  c  F.  =  2t  (t-1)(t-3) 

11  *  12  K 

(6.6) 

c  E.  +  c  F.  =  (b/a)  t2  (t-1)2  (t-3) 

21  K  22  K 

(6.7) 

where  the  c. .  are  given  in  Table  2  and  are  functions  of  the  ratio  of  coning 

*  vi 

frequency  to  spin  (t),  reduced  fineness  ratio  (f*  =  c/ka),  and  fill  ratio  f. 


Table  2.  Coefficients  in  the  Equations  (6.6,  6.7) 
That  Determine  Ek  and  Fk  for  §a  =  e  =  0 


ii 


12 


2  1 


2  2 


X2  = 


[(t  -  1)  X/f*]0x  * (X/f*)  -  E01(X/f*) 

[(t-1)  X/f*]Y  ‘(i/f*)  -  2Y  (X/f*) 

(t-1 )(x/f*)(b/a)  J  1  (xb/f*a) 
i 

+  '( x 2 - 1 ) (t-3)-2  J  (Xb/f*a) 

L  J  i 

(x-1 )(Xb/f*a)  Y  '  (Xb/f*a) 

i 

+  [(t2-1  ) (x-3) -2]  Yi  (xb/f*a) 

/  n  \  2  F  3  +  2t  -  t2 
\  2  )  1  -  2t  +  t2 


~The  k  =  0  mode  is  only  present for  liquid  payload  offset  h  f  0.  This  small 
moment  tern  is  included  in  the  computer  program.  Since  h  is  zero  in  all 
available  experimental  data }  the  zei'o  mode  brill  not  be  considered  further  in 
this  report . 
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For  fixed  reduced  fineness  ratio  and  fill  ratio,  the  determinant  of 
Eqs.  {6.6  -  6.7)  is  zero  for  certain  values  of  t  which  are  the  eigenvalues  of 
the  system,  xkn.  In  the  vicinity  of  these  eigenvalues,  the  parameters  Ek  and 

can  be  approximated  as  poles. 


E .  w 

k  '  T  •  Tkn 


(6.8) 


F  -  W 

k  ’  T  -  Tkn 


(6.9) 


Spiral  motions  are  represented  by  complex  values  of  s  while  pure  imagi¬ 
nary  values  of  s  correspond  to  constant  amplitude  coning  motion.  Eqs.  (6.8  - 
6.9)  can  now  be  extended  to  approximate  the  response  to  spiral  motions  in  the 
vicinity  of  the  eigenfrequencies  since  the  parameters  are  analytic  functions. 


E 


k 


Ek  (Tkn 
-It 


(6.10) 


1  Fk  <*kn> 


1  T 


kn 


(6.11) 


Stewartson  assumed  that  near  an  eigenfrequency  the  pressure  was  dominated  by 
that  mode  and  computed  the  total  pressure  moment  from  that  assumption. 


[f*  R  (f*.  f,  tkn)]2 
mp  2irk2  (s  -  i  xkn) 


(6.12) 


Stewartson3  and  later  authors16’17  have  constructed  tables  of  xkn  and  R  as 
functions  of  f*  and  f. 


A  portion  of  one  of  these  tables  for  a  fully  filled  shell  is  given  in 
Table  3.  As  can  be  seen  from  the  table,  R  decreases  rapidly  with  increasing 
n.  Since  the  moment  varies  as  R2,  eigenvalues  for  n  >  3  are  of  little  impor¬ 
tance  for  estimating  the  liquid  moment.  Eq.  (6.12)  also  shows  the  decay  of 
the  moment  with  increasing  k  values.  For  these  reasons,  our  computer  code 
only  considers  the  first  ten  k  modes  (k  =  1,  3,  5,  ...,  19). 


16.  B.  G.  Karpov ,  "Dynamics  of  Liquid-Filled  Shell.  Aids  for  Designers: 
(a)  Milne’s  Graphs  ( b )  Stewai'tson  ' a  Tables,"  BRL  Memorandum  Report 
1477 ,  May  196S.  AD  410484. 

17.  J.  T.  Fmsier,  " Dynamics  of  Liquid-Filled  Shell:  Aids  for  Designers 
BRL  Memoixmdm  Report  1892,  December  1967.  AD  665356. 
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Table  3.  f*  and  R  for  100%  Filled  Cylinder  (f  =  1) 


n 

=  1 

n 

=  2 

n 

=  3 

n 

=  4 

Tkn 

f* 

R 

f* 

R 

f* 

R 

__  .  . 

f* 

R 

0 

.9949 

0 

.4780 

0 

.3103 

0 

.2291 

0 

.01 

1  .0064 

.0144 

.4842 

.00172 

.3144 

.00048 

.2322 

.00019 

.02 

1.0180 

.0290 

.4904 

.00349 

.3186 

.00097 

.2353 

.00039 

.03 

1.0297 

.0440 

.4968 

.00532 

.3228 

.00148 

.2384 

.00060 

.04 

1.0417 

.0592 

.5032 

.00721 

.3271 

.00201 

.2416 

.00081 

.05 

1.0538 

.0747 

.5098 

.00915 

.3315 

.00256 

.2449 

.00104 

.06 

1.0661 

.0905 

.5165 

.01116 

.3360 

.00312 

.2482 

.00127 

.07 

1.0786 

.1065 

.5233 

.01323 

.3405 

.00370 

.2516 

.00150 

.08 

1.0913 

.1229 

.5302 

.01536 

.3451 

.00431 

.2550 

.00175 

.09 

1.1042 

.1395 

.5372 

.01755 

.3498 

.00493 

.2585 

.00200 

.10 

1.1173 

.1564 

.5443 

.01981 

.3546 

.00558 

.2621 

.00227 

.11 

1.1306 

.1736 

.5516 

.02214 

.3595 

.00624 

.2657 

.00254 

.12 

1.1442 

.1911 

.5591 

.02453 

.3645 

.00693 

.2695 

.00282 

.13 

1.1579 

.2088 

.5666 

.02700 

.3695 

.00764 

.2732 

.00311 

.14 

1.1720 

.2269 

.5743 

.02954 

.3747 

.00837 

.2771 

.00341 

.15 

1.1862 

.2452 

.5822 

.03215 

.3800 

.00913 

.2810 

.00372 

.16 

1.2008 

.2638 

.5902 

.03484 

.3854 

.00992 

.2851 

.00405 

.17 

1.2156 

.2826 

.5983 

.03761 

.3908 

.01073 

.2892 

.00438 

.18 

1  .2307 

.3018 

.6067 

.04046 

.3965 

.01156 

.2934 

.00472 

.19 

1.2460 

.3213 

.6152 

.04339 

.4022 

.01242 

.2977 

.00508 

.20 

1.2617 

.3410 

.6239 

.04640 

.4080 

.01332 

.3020 

.00545 

.21 

1.2777 

.3610 

.6327 

.04951 

.4140  ■ 

.01424 

.3065 

.00583 

.22 

1.2940 

.3813 

.6418 

.05270 

.4201 

.01519 

.3111 

.00623 

.23 

1.3107 

.4019 

.6511 

.05599 

.4264 

.01617 

.3158 

.00664 

.24 

1.3277 

.4227 

.6605 

.05937 

.4328 

.01719 

.3206 

.00706 

.25 

1.3450 

.4439 

.6702 

.06284 

.4394 

.01824 

.3255 

.00750 

.26 

1.3628 

.4653 

.6801 

.06642 

.4461 

.01933 

.3305 

.00795 

.27 

1.3809 

.4871 

.6903 

.07010 

.4529 

.02045 

.3357 

.00842 

.28 

1.3995 

.5091 

.7007 

.07389 

.4600 

.02161 

.3410 

.00891 

.29 

1.4185 

.5314 

.7113 

.07779 

.4672  j 

.02281 

.3464 

.00941 

.30 

1  .4379 

.5540 

.7222 

.08180 

■ 

.4746  | 

.02406 

.3520 

.00994 

The  major  difficulty  with  the  inviscid  Stewartson  eigenfrequency  theory 
is  that  it  predicts  much  too  large  liquid  moments  near  the  eigenfrequencies. 
Wedemeyer  introduced  a  viscous  boundary  layer  theory  and  sought  to  predict  the 
resulting  liquid  moment  by  a  simple  manipulation  of  Stewartson's  tables5.  Bv 
a  very  clever  approach,  he  showed  that  viscous  eigenvalues,  s^n»  could  be 

computed  from  Stewartson's  table  of  inviscid  eigenfrequencies,  x<,  *  by  the 

relations  kn0 


kn  '  (ckn  *  ()  Tkn 

(6.13) 

kn  Tkno  +  iTkn 

(6.14) 
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wh  i  re 


ekn  Tkn  +  iATkn  =  1 


9t 


kn 


9t 


3f*  ^f*6a  “  k  6c)  +  2 


kn 


9f 


(f  -  1)  S3 


He  then  replaced  i in  the  denominator  of  Eq.  (6.12)  wi  „h  the  s^n  of  Eq. 

(6.13)  to  obtain  an  excellent  approximation  of  the  pressu.  _>  moment  for  large 
Reynolds  numbers. 


The  maximum  value  of  the  pressure  side  moment  occurs  for  t  near  the 
Stewartson  eigenfrequency  and  is  approximately 


-  (f*R)2 _ 

2,k2  V  [ ekn  -  £] 


(6.15) 


For  e  =  0  and  constant  f*  the  maximum  SW  side  moment  coefficient  varies 
inversely  with  k 2 and,  hence,  it  varies  as  Re  -.  Its  dependence  on  k  is 

somewhat  more  complicated  since  k  5  varies  as  k_i  for  constant  f*.  The 

maximum  side  moment  coefficient  should  vary  with  k  by  a  factor  between  k 
and  unity. 

Equations  (6.3  -  6.5)  have  been  coded  for  N  =  19  by  Bradley18  and  a  num¬ 
ber  of  computations  made  for  a  variety  of  values  of  c/a,  f,  and  Re.  To  il¬ 
lustrate  his  results,  a  series  of  calculations  have  been  made  for  near 
.07,  e  =  0,  .02,  and  Re  =  500,000  and  15,000.  Table  3  shows  that  a  suitable 
value  of  f*  to  obtain  Ty  near  .07  is  1.08.  For  the  first  three  k-modes,  the 
corresponding  fineness  ratios,  c/a,  are  1.08,  3.24  and  5.40. 


Figures  1  and  2  show  the  side  moment  for  fineness  ratio  1.08  at  two 
Reynolds  numbers.  The  ratio  of  the  square  roots  of  the  two  Reynolds  numbers 
is  5.8  while  the  ratio  of  peak  CLSM's  for  c  =  0  at  the  two  Reynolds  numbers 

is  5.6.  Notice  the  strong  sensitivity  of  the  side  moment  to  damping  per 
cycle  for  the  higher  Reynolds  number.  This  sensitivity  is  considerably  re¬ 
duced  for  Rg  of  15,000. 

Figures  3-4  show  similar  curves  for  a  fineness  ratio  of  3.24.  The  max¬ 
imum  side  moment  for  the  higher  Reynolds  number  is  reduced  by  a  factor  of  5 
from  that  of  Figure  1.  Since  k  is  3,  the  predicted  range  of  this  factor  is 
1  to  9.  The  sensitivity  to  damping  per  cycle  is  quite  similar  to  that  shown 
by  Figures  1-2. 


18,  J,  V.  Bradley,  "Calaulations  of  liquid  Pay  lead  MceiMt , 
Report  in  prepatMtion. 


Mcrtorandur. 
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Finally,  Figures  5-6  show  the  side  moment  for  a  fineness  ratio  of  5.40. 
Here  the  reduction  of  peak  side  moment  in  Figure  5  relative  to  Figure  1  is  by 
a  factor  of  12,  which  also  lies  in  the  predicted  range  of  1  to  25.  The  de¬ 
pendence  on  damping  per  cycle  and  Reynolds  number  is  quite  similar  to  that  of 
the  preceding  figures. 

Since  the  in-plane  moment  only  affects  frequency,  it  is  of  much  less 
interest  than  the  side  moment.  It  is,  of  course,  available  from  Eqs.  (6.3  - 
6.5).  Figure  7  is  an  example  of  the  in-plane  moment  coefficient  for  a  Rey¬ 
nolds  number  of  500,000  and  fineness  ratio  of  3.24.  For  zero  x,  it  is  quite 
near  the  frozen  value  of  i,  undergoes  a  disturbance  near  the  eigenfrequency, 
and  asymptotically  approaches  zero. 

The  side  moment  has  contributions  from  the  two  flat  endwalls  as  well  as 
from  the  cylindrical  lateral  wall.  The  ratio  of  the  two  contributions  is 
shown  in  Figures  8-9  for  the  two  Reynolds  numbers  of  500,000  and  15,000.  We 
see  that  these  contributions  are  opposing  and  roughly  equal.  Indeed  for  a 
fineness  ratio  of  1.08  (k  =  1)  the  lateral  contribution  is  only  20%  larger  in 
magnitude  than  the  endwall  contribution!  Thus  the  side  moment  is  the  dif¬ 
ference  of  two  nearly  equal  quantities,  and  a  small  change  on  one  wall  could 
have  a  large  effect  on  the  side  moment. 


7.  WALL  SHEAR  MOMENT 


In  addition  to  the  pitch  and  yaw  moment  due  to  pressure  on  the  walls  of 
the  liquid  container,  moments  due  to  viscous  wall  shear  are  present.  These 
can  be  computed  from  the  derivatives  of  the  viscous  velocity  components  of 
Section  4.  The  liquid  moment  coefficient  due  to  shear  on  the  cylindrical  lat¬ 
eral  wall  is 


where 


mv£  =  (2ttK  Re)'1  e"s* 
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The  velocity  derivatives  in  Eq.  (7.2)  can  be  computed  from  Eqs.  (4.9  -  4.10) 
and  the  resulting  equation  has  been  coded  by  Bradley.16 

The  wall  shear  moment  coefficient  on  the  forward  flat  end  wall  is 


mvel  =  Re^_1  e"^  (ac^"1/*a  f  e18m*velr  de  dr 


where 


(7.3) 


m*vel  "  (h  +  c) 


» '1e 
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A  similar  expression  applies  for  the  rearward  endwall .  The  sum  of  these 
moment  coefficients  has  the  simple  form: 

A 

X  =  1 


mve  =  (2aK  Re) 


/ar 


-  iv  ) 
sv  sv' 


x  *  -1 


rdr  +  (h/c)2  mveh  (7.4) 


where 


mveh  *  <2K  Re> 
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The  velocity  derivatives  in  Equation  (7.4)  can  be  found  by  use  of  Equation 
(4.17)  and  the  results  have  also  been  coded  by  Bradley.18 

Since  all  the  velocity  derivatives  are  proportional  to  1  and  is 

i  a  a 

proportional  to  Re'5,  the  viscous  liquid  moment  coefficient  itself  varies  as 

Re-^  and  is  important  for  low  Reynolds  number.  In  Figure  10  the  wall  shear 
side  moment  coefficient  is  given  for  Re  =  15,000  and  our  three  sample  fine¬ 
ness  ratios.  Comparing  this  component  with  the  pressure-induced  liquid  side 
moment  coefficients  of  Figures  2,4,  and  6,  we  see  that  the  maximum  wall 
shear  component  is  from  10%  to  35%  of  the  wall  pressure  component. 


In  Figures  11-13,  the  total  side  moment  coefficients  for  Re  =  1,000  and 
e  =  0  are  compared  with  their  pressure  components  for  c/a  =  1.08,  3.24,  and 
5.40.  The  differences  between  these  curves  are  the  wall  shear  components  and 
we  see  that  these  components  are  quite  important  and  must  be  computed  for  low 
Reynolds  flows. 
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D'Amico  and  Miller  show  a  very  interesting  effect  of  very  low  Reynolds 
numbers19.  In  a  special  gyroscope  experiment20,  Miller  forced  a  spinning 
cylinder  filled  with  liquid  to  precess  at  an  angle  of  20°  and  measured  the 
despin  moment.  The  fineness  ratio  was  4.29,  t  varied  between  .12  and  .25. 
and  liquids  with  kinematic  viscosities  varying  between  1  and  106  centi stokes 
were  tested.  As  can  be  seen  from  Figure  14,  the  despin  moment  varied  from  a 
small  value  for  water  (v  =  1  cs)  to  a  value  thirty  times  bigger  for  corn 
syrup  (v  =  2  x  10s  cs).  The  Reynolds  number  of  this  peak  was  about  10. 

D'Amico  and  Miller  conjectured  that  this  large  despin  moment  would  be 
associated  with  a  large  side  moment  which  would  produce  flight  instabilities. 
Flight  tests  were  made  and  this  otherwise  unexpected  instability  was  observed. 
We  now  can  use  the  theory  of  this  report  to  estimate  the  liquid  side  moment 
for  Miller's  cylinder  and  Reynolds  number  as  low  as  100,  which  is  probably 
the  extreme  lower  bound  of  validity  of  a  boundary  layer  theory. 


In  Figure  15a,  is  plotted  versus  t  for  Miller's  cylinder  with 

Re  =  10s.  The  four  local  maxima  on  these  curves  are  caused  by  four  eigen- 
frequencies,  Tnk.  By  use  of  Table  3  these  eigenfrequencies  can  be  identified. 

In  order  of  increasing  t  their  (k,  n)  mode  numbers  are  (15,  4),  (11,  3), 

(7,  2)  and  (13,  4).  As  we  would  expect,  the  largest  maximum  has  n  =  2  while 
the  two  quite  small  maxima  have  n  =  4. 


In  Figures  I5b-~15e,  C^SM  is  computed  for  Re  =  TO5,  104,  1 03 ,  102.  The 

first  effect  of  decreasing  Reynolds  number  is  to  decrease  the  size  of  the 
maxima  associated  with  the  eigenfrequencies.  Next  we  see  that  the  average 
level  of  the  side  moment  coefficient  curves  increases  with  decreasing  Rey¬ 
nolds  number.  In  Figure  16,  is  plotted  versus  Re  for  t  =  .10,  .15,  .20, 

.25.  To  facilitate  comparison  of  these  curves,  we  normalized  each  side 
moment  coefficient  by  its  values  at  Re  =  10s.  We  see  that  the  side  moments 
increase  to  maximum  values  18  -  33  times  their  values  for  wate**  and  these 
maxima  occur  around  Re  =  300.  This  striking  qualitative  agreement  with  the 
D | Ami co-Mi lie:*  conjecture  is  very  exciting,  and  theoretical  work  on  pre¬ 
dicting  the  despin  moment  at  low  Reynolds  number  is  being  given  a  much 
greater  emphasis  as  a  result. 


1 9.  X.  F.  D'Amico  and  M.  C.  Mil  lev,  "Flight  Instability  Produced  by  c 
Rapidly  Spinning,  Highly  Viscous  Liquid, "  Journal  of  Sj>accci\ift  and 
Rockets,  Vol .  16,  January -February  1979,  pp.  62-64, 

20,  M.  C .  Miller,  "Flight  Instabilities  of  Spinning  Projectiles  Raving  Ron- 
rigid  Payloads,"  Journal  of  Guidartce,  Control ,  and  Dynamics,  Vol.  5, 

pp .  151-157,  Mai'ch-April  1962. 
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8.  EXPERIMENTAL  RESULTS 


In  1980  the  first  wall  pressure  measurements  in  a  precessing  and  spin¬ 
ning  liquid-filled  cylinder  were  made  by  Whiting21.  For  the  100%  filled  cyl¬ 
inder,  he  compared  his  measurements  with  theoretical  calculations  by  Gerber 
et  al22.  Gerber's  theory  was  developed  for  fully-filled  cylinders  only  and 
was  used  to  compute  the  viscous  influence  of  the  lateral  wall  exactly  without 
the  use  of  a  boundary  layer  approximation.  Its  results  for  fully-filled  cyl¬ 
inders  should,  therefore,  be  better  than  that  of  this  report  when  they  differ. 
For  the  Reynolds  numbers  of  the  Whiting  tests,  they  did  not  differ  signifi¬ 
cantly;  and  thus,  the  good  experimental  agreement  Whiting  got  for  Gerber's 
calculations  also  applies  to  our  theory.  In  two  cases,  however,  Whiting 
measured  the  pressure  on  the  flat  end  wall  of  a  partially-filled  cylinder 
(f  =  .92).  Comparisons  of  these  measurements  with  the  theoretical  prediction 
of  Equation  (6.3)  is  given  in  Figures  17a-17b.  The  agreement  with  theory  is 
quite  satisfactory. 

In  most  gyroscope  experiments6’23  the  yaw  growth  rates  and  coning  rates 
are  measured  for  a  variety  of  test  conditions.  In  all  experiments  the  center 
of  mass  is  located  at  the  pivot  point  so  that  the  gyroscopic  stability  factor 
is  infinitely  large  and  Equations  (2.16,  2.18)  for  frequency  and  damping  be¬ 
come: 


=  0  [l  +  (m^/g  C 

*  (v!A) 


=  0 


-1 


(2x/o  -  1)  CLSM  (x,e) 


(8.1) 

(8.2) 


We  first  consider  D'Amico  and  Rogers'23  measurements  for  a  cylinder  with 

fineness  ratio  of  1.042  for  which  x  =  .040.  The  frequency  was  changed  by 

1 1 


Comparisons  of  theory  with  these  data  are  shown  in  Figures  18a-18b.  Agreement 
is  fair,  but  there  appears  to  be  a  systematic  bias.  Theoretical  curves  for 
different  fineness  ratios  in  steps  of  .001  were  computed  and  the  best  fits  are 
shown  as  dashed  curves  for  c/a  =  1.047  and  1.048  in  the  respective  figures. 

An  effective  fineness  ratio  .5%  greater  than  the  measured  value  is  not  unrea¬ 
sonable  and  gives  excellent  experimental  agreement.  Figure  19  compares  the 
complete  side  moment  coefficient  with  pressure  side  moment  coefficient  and  its 
Stewartson-Wedemeyer  approximation  for  c/a  =  1,048  and  Re  =  2400,  and  we  see 
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that  the  complete  moment  coefficient  which  gave  such  excellent  experimental 
agreement  is  quite  different  from  the  other  more  approximate  curves. 

The  remaining  three  sets  of  experiments  are  summarized  by  Whiting  and 
Gerber6,  and  all  involve  very  similar  fineness  ratios;  i.e.,  c/a  =  3.149, 
3.013,  3.077.  The  Stewartson  eigenfrequency  for  the  first  set,  which  was  the 
only  one  that  had  100%  filled  cylinders,  is  x^  "  .047.  The  appropriate 

eigenfrequency  for  the  other  two  is  x^,  but  the  proper  values  cannot  be  ob¬ 
tained  from  Table  3. 


In  the  first  two  sets  of  experimental  data  presented  oy  Whiting  and 

Gerber,  the  frequency  was  changed  by  varying  I  and  I  .  For  this  procedure, 

x  y 

Equation  (8.2)  gives  e  as  a  function  of  x  and  I  and  thus  a  simple  theoretical 

curve  cannot  be  plotted.  What  was  done  was  to  compute  theoretical  values  for 
each  experimental  pair  of  values  (x,  Ix),  plot  these  values  as  points  in 

Figures  20a-b  and  21a-d  and  connect  the  points  with  straight  Vnes.  An  exam¬ 
ination  of  these  figures  shows  fair  to  good  agreement  between  theory  and 
experiment.  It  is  interesting  to  note  that  agreement  with  experiment  for 
Re  =  520,000  (Fig.  20a)  can  be  considerably  improved  by  a  0.1%  change  in  fine¬ 
ness  ratio.  The  side  moment  coefficients  for  the  lowest  Reynolds  number  of 
each  set  (namely,  9000  and  5200)  are  given  in  Figures  22-23. 

In  the  final  experiment  to  be  considered,  the  fill  ratio,  f,  was  varied 
and  the  frequencies  and  yaw  damping  rates  were  measured.  Here,  too,  the 
theoretical  yaw  damping  rate  is  a  function  of  two  variables  -  fill  ratio  and 
frequency  -  and  must  be  represented  by  individually  computed  points  connected 
by  line  segments.  Although  the  agreement  for  Re  =  520,000  given  in  Figure 
24a  is  quite  good,  the  situation  for  Re  =  5,200  in  Figure  24b  is  poor. 

9.  CENTRAL  ROD 


In  1969  Frasier7’6  extended  the  SW  theory  to  the  fully-filled  cylinder 
with  a  central  rod.  This  had  the  effect  of  replacing  the  free  surface  bound¬ 
ary  condition  of  Eq.  (5,8)  with  an  inner  lateral  surface  condition  of  the 
same  form  as  the  outer  lateral  surface  (Eqs.  (4.9  -  11)).  For  a  rod  with 
radius  d,  these  conditions  are 
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Equation  (9.3)  can  be  used  to  derive  an  inner  boundary  condition  to  substi¬ 
tute  for  Eq.  (5.21). 


21  (1  -  a6a/d)  Rk(d)  -  [s  -  i  (1  +  2a6a/d)J  d  R * k(d) 


-  (s  -  1)  6a  adR"(d)  =  2(d/a)  aks  (s  -  i)(s  -  3i)  (9.4) 

Using  boundary  conditions  (5.20)  and  (9.4)  for  6a  *  0,  Frasier  and  Scott7 

a 

calculated  tables  of  Inviscid  elgenfrequencles,  xnk. 

These  frequencies  are  functions  of  the  reduced  fineness  ratio  and  rod- 
ded  fill  ratio  f^. 


Tkn  =  Tkn  (f*’  V 


(9.5) 


where 


fd  =  1  -  d*/a2 


The  revised  Wedemeyer  relation  for  the  viscous  eigenvalue  is  simply 
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Table  4  is  a  sample  table  of  tkn,  f*,  and  R  for  f^  *  .98  (d/a  =  .14).  A 
similar  table  of  tkr),  f*,  and  R  for  fill  ratio  (f)  of  .98  is  given  as  Table  5 
fAr  comparison  purposes. 

Two  more  moment  coefficient  terms  must  be  added  to  E^.  (6.1).  These 
are  due  to  the  pressure  on  the  rod,  m  ,  and  the  wall  shear  on  the  rod,  myr. 
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Table  4.  f*  and  R  for  Cylinder  With  Rod  (f^  =  .98) 


n 

=  1 

n 

=  2 

n 

=  3 

n 

=  4 

Tkn 

f* 

R 

f* 

R 

f* 

R 

f* 

R 

.00 

.9028 

.0000 

.4102 

.00000 

.2615 

.00000 

.1921 

.00000 

.01 

.9141 

.0127 

.4160 

.00092 

.2652 

.00039 

.1948 

.00008 

.02 

.9256 

.0257 

.4219 

.00189 

.2690 

.00080 

.1975 

.00017 

.03 

.9373 

.0389 

.4279 

.00290 

.2728 

.00122 

.2003 

.00027 

.04 

.9491 

.0525 

.4339 

.00395 

.2766 

.00165 

.2031 

.00036 

.05 

.9611 

.0664 

.4401 

.00505 

.2806 

.00210 

.2059 

.00046 

.06 

.9733 

.0806 

.4464 

.00620 

.2846 

.00257 

.2089 

.00057 

.07 

.9856 

.0951 

.4528 

.00740 

.2887 

.00305 

.2118 

.00068 

.08 

.9981 

.1100 

.4593 

.00865 

.2929 

.00355 

.2149 

.00079 

.09 

1.0109 

.1251 

.4659 

.00995 

.2971 

.00407 

.2180 

.00091 

.10 

1.0238 

.1405 

.4726 

.01130 

.3015 

.00461 

.2211 

.00103 

.11 

1.0369 

.1562 

.4795 

.01271 

.3059 

.00516 

.2244 

.00116 

.12 

1.0503 

.1723 

.4865 

.01418 

.3104 

.00574 

.2276 

.00130 

.13 

1.0638 

.1886 

.4936 

.01571 

.3150 

.00633 

.2310 

.00144 

.14 

1.0776 

.2053 

.5009 

.01729 

.3197 

.00694 

.2344 

.00159 

.15 

1.0916 

.2223 

.5083 

.01894 

.3245 

.00758 

.2379 

.00174 

,16 

1.1059 

.2395 

.5159 

.02066 

.3294 

.00824 

.2415 

.00190 

Table  5.  f*  and  R  for  982  Filled  Cylinder  (f  = 

.98) 

n 

=  1 

n 

=  2 

n  = 

3 

n 

=  4 

Tkn 

f* 

R 

f* 

R 

f* 

R 

f* 

R 

.00 

.9942 

.0000 

.4753 

.00000 

.3053 

.00000 

.2224 

.00000 

.01 

1.0057 

.0144 

.4814 

.00171 

.3094 

.00047 

.2254 

.00018 

.02 

1.0173 

.0290 

.4877 

.00346 

.3135 

.00095 

.2284 

.00037 

.03 

1.0291 

.0439 

.4940 

.00527 

.3176 

.00145 

.2315 

.00057 

.04 

1.0410 

.0591 

.5005 

.00713 

.3219 

.00197 

.2346 

.00077 

.05 

1.0532 

.0746 

.  5070 

.00904 

.3262 

.00250 

.2378 

.00098 

.06 

1.0656 

.0904 

.5137 

.01100 

.3307 

.00306 

.2410 

.00119 

.07 

1.0781 

.1064 

.5205 

.01302  : 

.3352 

.00363  j 

.2443 

.00141 

.08 

1.0909 

.1227 

.5275 

.01510 

.3398 

.00423 

.2*77 

.00164 

.09 

1 . 1039 

.1393 

.5346 

A1  70A 

•  N/  «  i 

7  A  A 

AAJ  ft  * 

*  * 

.2511 

AA  *  AA 

.W1W 

.10 

1.1171 

.1562 

.5419 

.01944 

.3492 

.00548 

.2547 

.00212 

.11 

1.1306 

.1733 

.54  92 

.02170 

.3541 

.00614 

.2583 

.00237 

.12 

1.1443 

.1907 

.5567 

.02402 

.3591 

.00682 

.2619 

.00263 

.13 

1.1583 

.2084 

.5644 

.02641 

.3642 

.00752 

.2657 

.00290 

.14 

1.1725 

.2264 

.5722 

.02887 

.3694 

.00825 

.2695 

.00318 

.15 

1.1871 

.2447 

.5802 

.03139 

.3747 

.00901 

.2734 

•00347 

.16 

1.2019 

.2632 

.  5884 

.03399 

.3802 

.00979 

.2774 

-00376 
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Frasier  ran  gyroscopic  experiments  for  a  rodded  cylinder  with  c/a  = 
2.864,  =  .977  (d/a  =  .15)  and  three  Reynolds  numbers.  Figures  25a  -  25c 

show  his  data  and  the  improved  SW  theoretical  prediction.  The  agreement  is 
rather  good  although  the  peak  CLSfyj  is  not  predicted  very  well.  In  Figure  26 

the  side  moment  coefficient  for  frf  =  .977  is  compared  with  that  for  f  =  .977. 

The  side  moment  for  rodded  cavity  shows  an  eigenfrequency*  at  .059  while  that 
for  partially-filled  cavity  shows  no  eigenfrequencies  and  is  a  small  negative 
value.  Thus,  we  s  1  that  the  presence  of  an  inner  cylindrical  wall  can  have 
a  very  strong  effect  on  the  side  moment.  The  insertion  of  a  cylindrical 
partition  to  improve  the  stability  of  liquid-filled  shell  has  been  proposed 
by  Frasier  and  D'Amico5,  and  D'Amico24  has  experimentally  studied  the  side 
moment  during  transition  from  a  free  surface  to  a  fully-wetted  central  rod. 


24,  W.  P.  D'Amico ,  "Dynamics  of  Liquid  Filled  Shell:  Liquid-Central  Burster 
'  Interference 3"  BEL  Memorandum  Report  1985,  June  1969.  AD  855134. 

*  According  to  Table  is  .045  for  =  ,980  and  we  would  identify  this 

peak  at  .059  in  Fig.  26  to  be  caused  by  x 
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10.  OTHER  THEORETICAL  EXTENSIONS 


The  extended  SW  theory  of  this  report  has  a  number  of  restrictions. 

Some  of  these  restrictions  have  been  relaxed  by  the  work  of  a  number  of  re¬ 
searchers.  These  restrictions  include  those  requiring:  (1)  a  single  liquid 
filler:,  (2)  a  centrally  located  circular  cylinder  container;  (3)  small 
amplitude  motion;  (4)  steady-state  motion;  and  (5)  fully  spun-up  liquid. 

In  this  section  we  will  survey  the  work  in  these  five  areas. 


10.1  Single  Liquid  Filler 

In  1972  Scott25  considered  the  eigenfrequencies  for  an  inviscid  two- 
component  liquid.  He  included  the  possibility  of  partial  fill,  computed 
eigenfrequencies,  and  got  fair  experimental  agreement. 


10.2  Centrally  Located  Circular  Cylinder  Container 

Scott26  also  considered  the  inviscid  eigenfrequencies  and  moments  for 
eccentrically  located  fully  filled  circular  cylinders  and  showed  that  these 
frequencies  and  moments  are  the  same  as  those  for  a  centrally  located  cylinder. 
In  an  earlier  work  Wedemeyer27  derived  an  approximate  relation  for  inviscid 
eigenfrequencies  of  a  slightly  noncylindrical  cavity.  He  showed  that  the 
Stewartson  eigenfrequency  tables  could  be  used  through  the  use  of  an  average 
fineness  ratio. 


(c/a) 


av 


(10.1) 


Karpov28  made  a  number  of  gyroscope  experiments  that  showed  good  results  for 
this  Wedemeyer  concept. 


25.  W.  E.  Scott,  "The  Inerti.il  Wave  Frequency  Spectrun:  in  a  Cylindrically 
Confined,  Inviscid,  Incompressible ,  Two  Component  Liquid, "  BEL  Report 
1609,  September  1972.  AD  752459. 

26.  W  B.  Scott,  "The  Dynamical  Effect  of  Inertial  'Menu-van  the  Free  Flight 
Motion  of  a  Body  Containing  Several  Eccentrically  Located,  Liquid- 
Filled  Cylinders ,"  BRL  Report  1551,  September  1971.  AD  755565. 

27.  E.  H.  Wedemeyer,  " Dynamics  of  Liquid  Filled  Shell:  Ron-Cylindrical 
Cavity,"  BRL  Report  1526,  August-  1966.  AD  489899. 

28.  B.  G.  Karpov,  "Dyrur-irs  of  Liquid-Filled  Shell:  Resonance  in  Modified 
Cylindj'ical  Cavities,"  BRL  Report  1552,  August  1966.  AD  804828. 
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10.3  Small  Amplitude  Moticn 

Gyroscope  experiments  by  Scott  and  D'Amico29  have  shown  that  the  yaw 
growth  rate  changes  from  the  linear  SW  values  at  coning  angles  as  low  as  1°. 
Indeed,  pressure  coefficient  measurements  by  Whiting21  have  shown  nonlineari 
ties  for  coning  angles  as  low  as  .05°!!  The  Scott-D'Amico  data  also  showed 
a  shift  in  eigenfrequency  for  coning  angles  in  excess  of  1°.  Scott30  later 
derived  a  modified  fineness  ratio  that  was  fairly  good  in  predicting  this 
frequency  shift. 


10.4  Steady-State  Motion 

D'Amico  et  al14  made  liquid  pressure  coefficient  measurements  on  the 
endwall  of  a  spinning  cylinder  whose  coning  motion  was  impulsively  started 
and  found  that  cone-up  time  for  the  pressure  coefficient  to  reach  the  steady- 
state  SW  value  could  be  as  large  as  five  seconds.  He  made  an  estimate  of 
this  time  from  the  real  part  of  the  appropriate  and  found  good  agreement. 

Since  these  cone-up  times  are  a  significant  part  of  a  projectile's  flight 
time,  work  in  this  important  area  is  continuing. 


10.5  Fully  Spun-Up  Liquid 


This  area  has  received  much  more  attention  than  that  given  to  the 
preceding  four  areas.  Wedemeyer31  developed  a  very  simple  model  of  the 
spin-up  process  which  was  extended  by  Kitchens  et  al3  ‘3  .  Spin-up  times 
greater  than  the  cone-up  times  have  been  predicted  and  measured..  Karpov34 
made  use  of  Wedemeyer ‘s  suggestion  to  obtain  an  estimate  of  the  effect  of 
spin-up  on  the  liquid  side  moment.  This  very  approximate  result  is  new  being 


29.  W.  E .  Scott  and  W.  P,  D'Amico ,  "Amplitude-Dependent  Behavior  of  a  Liquid- 
Filled  Gyroscope,"  Journal  of  Fluid  Mechanics ,  Vol,  SO,  Part  4 ,  2972, 
pp.  751-758. 


30.  W.  E.  Scott,  " The  Large  Amplitude  Motion  of  a  Liquid-Filled  Gyroscope  and 
the  Non-Interaction  of  Inertial  and  Rossby  Waves,"  Journal  of  Fluid 
Mechanics ,  Vol.  72,  Pai't  4,  1975,  pp.  649-860. 

21.  E.  H.  Wedemeyer,  " The  Unsteady  Flo’J  Within  a  Spinning  Cylinder,''  Journal 
of  Fluid  Mechanics ,  Vol.  SO,  Part  2,  1964,  pp.  383-399.  (Sea  also  BFJ. 
Report  1.225,  October  1965 ,  Ap  431946.  • 


C.  W.  Kitchens,  Jr.,  "Fkman  Compatibility  Conditions  in  Wedemeyer  Spin-Up 
Model, "  Physics  of  Fluids,  Vol.  23,  Part  5,  Wav  7960,  pp.  1062-1064. 


33.  C.  W.  Kitchens,  Jr. ,  N.  Gerber,  and  R.  Sedney,  "Sigin  Decay  of  Lit juui- 
Filled  th'ojecti les,"  Journal  of  Spacecraft  and  Pcckpts,  Vol.  18, 
Noverber-December  1878,  pp.  348-254.  (See  atso~BRl  Report  1998,  Jute 
1977,  AD  A043275  and  BRL  Report  2026,  October  1977,  AD  A050311.: 

34.  D.  G.  Karpov,  " Dynamics  of  Liquid-Filled  Shell:  Instability  During  Spin- 
Up,"  B8L  Memorandum  Report  1629,  January  1965.  AD  463826. 
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replaced  by  the  current  efforts  of  Sedney  et  al 13,22  to  develop  a  very  re¬ 
fined  perturbation  analysis  for  computing  the  liquid  side  moment  during 
spin-up.  The  only  direct  pressure  measurements  during  spin-up  have  been  made 
by  Aldridge35’36,  and  these  were  for  the  simple  case  of  axi symmetric  oscilla¬ 
tions  and  not  the  three-dimensional  oscillation  induced  by  coning  motion. 


11.  SUMMARY 

1.  A  general  definition  of  the  1 i quid  moment  has  been  developed,  and  the 
expression  for  frequencies  and  damping  of  projectile  angular  motion  has  been 
obtained. 

2.  An  exact  pressure  moment  has  been  computed  for  the  Stewartson-Wedemeyer 
theory. 

3.  Wall  shear  effects  have  been  added  to  the  improved  SW  pressure  moment. 

4.  The  improved  theory  shows  a  decrease  in  the  size  of  eigenfrequency- 
associated  peaks  in  the  side  moment  with  decreasing  Reynolds  number. 

5.  The  average  level  of  the  side  moment,  however,  grows  with  decreasing 
Reynolds  number  to  a  peak,  in  good  qualitative  agreement  with  the  D* Ami  co- 
Miller  conjecture. 

6.  Good  agreement  with  all  available  published  experimental  data  has  been 
shown. 
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Pressure  Component  of  CLS|vj  vs  x  for  Re  =  500,000,  c/a  =  1.08, 
1 ,  e  =  0  and  .02 


Figure  2.  The  Pressure  Component  of  C.  CM  vs  x  for  Re  =  15,000,  c/a 
f  -  1,  e  =  0  and  .02  Lirq 


Figure  4.  The  Pressure  Component  of  C.  vs  r  for  Re  =  15,000,  c/a 
f  =  1,  e  =  0  and  .02  L;>" 
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Figure  6.  Ti  e  Ratio  of  the  Lateral  Pressure  Component  of  C.  <.,,1  to  the  EndwalT 
Pressure  Component  vs  t  for  Re  =  b00,000,  f  =  1 ,  c  -  0,  c/a  =  1.08 
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Figure  15b.  C,  <-M  vs  x  for  Re  =  10s.  e  =  0,  and  the  D' Ami co-Mi ller  Parameters: 


artially-Filled  Cavity  (f  =  .920),  Predicted  Curve 

a  (see  Ref  21,  Fig  lOg),  Re  =  80,000,  c/a  =  3.148, 
/a  =  .668 


Figure  18a.  e  vs  r  for  a  Fully-Fined  Cavity,  Predicted  Curve  and  D'Amico 
Data  (see  Ref  23,  Fig  5),  Re  =  12,400,  c/a  =  1.042  (Nominal) 
and  1.047  (Fitted),  m,  a2/Iv  =  .0833 


0.02 
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e  vs  t  tor  a  ruily-MJled  Cavity,  Predicted  Values  and  DAmico- 
Kitchens  Data  (Circles)  (see  Ref  6,  Fig  8b),  Re  =  9000, 
c/a  =  3.148,  m^a2/!  Variable 


Figure  21a.  e  vs  t  for  a  Partially-Fi  11  ed  Cavity  (f  =  .790),  Predicted  Values 
and  D'Amico  Data  (Circles)  (see  Ref  6,  Fig  6a),  Re  =  520,000, 
c/a  =  3.013,  m.  a2/I  Variable 


e  vs  t  for  a  Partially-Fined  Cavity  (f  =  .790),  Predicted 
and  D'Amico  Data  (Circles)  (see  Ref  6,  Fig  6b),  Re  =  40,00) 
c/a  =  3.013,  m.a2/I  Variable 


Figure  21d.  e  vs  r  for  a  Partially-Fi 
and  D'Amico  Data  (Circles 
c/a  =  3.013,  m,  az/I  Vari. 


Has 


Figure  24b.  z  vs  Fill  Ratio  f.  Predicted  Values  and  Karpov  Data  (Circles) 
(see  Ref  6,  Fig  4),  Re  =  5200,  c/a  =  3.077 


e  vs  t  for  a  Cavity  With  Rod  (f .  =  .977),  Predicted  Curve 
and  rosier  Data  (see  Ref  6,  Fi^  9a),  Re  =  520,000, 
c/a  =  2.866,  m,  a2/I  =  .0240 
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for  a  Cavity  Kith  Rod  ' f  .  =  .977),  Predicted  Curve  and 
Data  (see  Ref  6.  Fig  Qbj,  Re  =  173,000,  c/a  =  2.866, 

=  .0265 
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APPENDIX  A 


DERIVATIONS  OF  EQUATIONS  (3. 7-3. 9) 


For  single  mode  coning  motion,  the  earth-fixed  components  of  unit  vectors 
along  the  non-rotating  aeroballistic  axes  have  the  following  simple  form 
accurate  to  the  first  order  in  K-. 


ej  =  (1,  -Kj  cos  -Kj  sin  $j)  (Al) 

et  =  (Ki  cos  1,  0)  (A2) 

y  j  j 

e^  =  (Kj  sin  4>j»  0,  1)  (A3) 


The  vector  from  the  projectile's  center  of  mass  to  any  point  on  the 

/V 

projectile  can  be  given  by  aeroballistic  cylindrical  coordinates  (x,  r,  e)  and 

can  then  be  related  to  earth-fixed  coordinates  (xe,  ye,  ze)  by  the  following 
equation: 


(xe»  ^e»  ze^ 


X  Vr 


cos  8  e  +  sin  0  e 

S  2 


(A4) 


If  we  now  introduce  earth-fixed  cylindrical  coordinates  (xe,  re,  0e),  the 
three  component  equations  of  vector  Equation  (A4)  are: 


xe  =  x  +  r  Kj  cos  Uj  -  0)  (A5) 
re  cos  ee  =  r  cos  8  -  x  Kj  cos  $j  (A6) 
re  sin  9e  =  r  sin  e  -  x  Kj  sin  ^  (A7) 


Since  the  earth-fixed  cylindrical  coordinates  will  be  used  throughout  this 
report,  and  missile-fixed  cylindrical  coordinates  are  never  used,  we  can  omit 
the  subscript  "e"  without  any  ambiguity  problem  and  will  do  so  as  a  con¬ 
venience.  Equation  (A6)  can  be  multiplied  by  sin  6,  Equation  (A7)  by  cos  0, 
and  the  results  subtracted  to  yield: 


A. 

sin  (o  -  0)  =  -(x  Kj/r)  sin  ($j  -  e) 


(AS) 
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or 


0  =  3  +  R  < i (x/r)  K  es*  -  16 


(A9) 


Next  we  square  Equations  (A6  -  A7)  and  add  the  results. 


r2  =  r2  -  2  x  r  Kj  cos  (<|>j  -  6)  +  Kj 


2X2 


(A10) 


or 


(All) 


Finally,  Equation  (A9)  can  be  used  to  obtain  a  revised  version  of  Equation 
(A5). 


x 


x  +  R  <  r  K  es$  "  i0 


(A12) 


For  any  fixed  point  on  the  projectile,  x  =  r  -  0,  e  =  Its  velocity 
in  earth-fixed  cylindrical  coordinates  can  be  computed  by  differentiating 
Equations  (A9,  All,  A12). 


«x=R  j  $  (s  -  i)rK  es$  ”  j 

(A13) 

°  r  a  -R  1 1  (s  -  1)xK  es$  ie  | 

(A14) 

=  r6  =  r$  +  R  |  i$  (s  -  i)xK  es$  “  i9| 

(A1S) 
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APPENDIX  B 


SOLUTION  OF  BOUNDARY  LAYER  EQUATIONS 

For  an  incompressible  fluid  with  constar.c  viscosity,  the  Navier-Stokes 
and  continuity  equations  in  cylindrical  coordinates  are: 


3V  ,,  3 V  Va  3 V  V2  3 V 
_ r  V  _ r  _0  _ r  _9  V  _ r 

3t  r  3r  r  36  "  r  x  3x 


r 

»L3r  L 


^  2  SV» 
?c  Vr"  r7  '  F  sT. 


!!etvr!\t\!!etVetvx!!e 

at  ar  r  36  r  3x 


r  y  ^y 

«  -  —  If  +  V  v2  V  -  -£  +  -L-  — r 
p  r  30  [_  c  6  r7  r2  30  J 


3 V  ,,  3 V  3 V  ..3V  .  , 

rr  ♦  Vr  r1  +  -  rr  ♦  Vx  ~  3  -  1  ^  ♦  v  V 


3r  r  36 


3x  p^3x  c  x 


where 


3V  V  ,  DV  3V 

^_n  +  j:  +  l_i  +  _j<c  o 

3r  r  r  DO  3x 


?  „  32  ,  1  3  A  32  I  D2 

c  Dr‘  r  Dr  Dx2  r*  301' 


Equations  (3.10  -  3.13)  can  be  substituted  in  Equations  (B 1  -  84),  and  prod¬ 
ucts  of  us>  v$,  w&  neglected  to  yield 


Dp  I"  a2v  2a 

(s  -  1)«s  -  2«s  **FS.  vRe-'|v-2  ,s  -  -Ts-  +  - 


*1  •  *1 
•  1W 

s 

rr_  J 
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iapt 

(s  -  i)u%  ♦  2vs  - 


yRe 


-l 


7o  "s 


ap, 

(s  '  i)us  +  a  ST  "  1rRe"1 5e  “s 


a(rv  )  au 

— - iw  +  r  =  0 

3r  5  sr 


a2wp  2a2iv„ 
s  -  s 

~rr'  r? 


where 


V0 


(B6) 

(B7) 

(B8) 


Next  we  assume  that  the  velocity  components  and  pressure  can  be  written  as  the 
sum  of  inviscid  and  viscous  terms.  The  inviscid  terms  satisfy  Equations  (B5  - 
B8)  for  Re-1  =  0  and  the  viscous  terms  satisfy  Equations  (B5  -  B8)  but  are  zero 
except  for  small  regions  near  the  walls  extending  a  distance  6  from  the  walls. 

We  will  make  the  usual  boundary  layer  assumptions  that  S^Re-1/2  and  derivatives 

normal  to  the  wall  vary  as  5_1  while  deviations  along  the  wall  are  of  order 
unity.  Positive  spin  will  be  assumed  (y=l )  since  the  effect  of  negat  ve  spin 
follows  from  Eq.  (2.8). 


=  u  ,  +  u 
si  sv 

(B9) 

V 

s 

=  V  .  +  V 

si  sv 

(BIO) 

w 

s 

=  w  .  +  w 
si  sv 

(Bll) 

ps 

=  p  .  +  n 

Hsi  Msv 

(B12) 

where  u  =  v  =w  =  p  =  0  far  from  wall  . 

sv  sv  sv  sv 

u$,  vs>  w&  must  satisfy  Eqs.  (3.7  -  3.9)  at  the  wa1!. 

•  •  usv 

e  (s  -  i)  (r/a)  K  -  . 

(B13) 

< 

<✓> 

< 

=  *(s  -  i)  (x/a)  K  -  vsi 

(814) 

w 

sv 

=  i(s  -  i)  (x/a)  K  -  w$1 

(B15) 
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Near  the  cylindrical  lateral  wall ,  r  “  a,  E,s.  (BS  -  B7)  become 


a2v 


a  =  2»sv  -  Is  '  ,)vsv  +  a2Re'1  ^ 

(s  -  i)w  =  a2Re-' 


sv 


a2"sv  -  2v  + 

7T  bV  r 


ar* 


a2u. 


(s-i)u  =  aZRe"1  ll|v  -  a 
^  ;  SV  3r2 


ap 


sv 

ax 


(B16) 

(B17) 

(B18) 


therefore,  v$v  is  sv  der  ^  and  can  be  neglected  in 

(B17).  (B16)  now  shows  that  psv  1S  °  sv 


Eq.  (Bi8) . 


a!!sv-2«sv 

ar 


(s  -  i)w$v  3  a2Re 


a2w 


-\ 


sv 


r2" 
a2u 


ar- 

i2, 


(s  *  ^usv  3  a2Re 


-i 


sv 


(B19) 


(B20) 


(B  21 ) 


ar* 


The  solutions  to  Eqs.  (B20  -  B21) 

^ (1  +  is)  (x/a)  K  ~  wsi] 

K 

1  ♦  i 


that  satisfy  Eqs.  (B13.  B15)  are: 
(r-a)/a6a 


u 


where 


sv 

4 

a 


-[(i  *  O  *  *  usi] 

1 

Re->/: 


(r-a)/a5. 


(B22) 

(B23) 


J 2(l+is) 


5ubstituting  Eqs.  (B22  -  B23)  in  -q 
near  the  lateral  «all  r  =  a. 

r  -  iva'l 

1(1  -  s)  t»/a)  H  -  VS1  •  a  1TJ 


(BB)  and  integrating,  we  can  obtain  v_v 

(r-a)/a6a  (82$) 


v  = 
sv  a 
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Eq.  (B24)  can  then  be  inserted  in  boundary  condition  (B14)  to  give  a  boundary 
condition  on  the  inviscid  radial  velocity  at  r  =  a. 


3V 


v„,  '  as 


si 


si  a  3r 


(1  -  s)  (x/a)  K 


(B25) 


Turning  now  to  the  endwalls,  x  =  (x-h)/c  =  +1,  similar  boundary  layer  size 
arguments  give  the  following  equations: 


(s  - 

i)vs.  -  2“sv 

=  a2Re_1  a# 

(826) 

(s  - 

,)wsv  +  2vsv 

32W 

■  a2«p-‘ 

(827 ) 

3psv 

3X 

a  0 

(B28) 

Next  Eq.  (B 26 )  is  multiplied  by  i  ano  both  added  and  subtracted  from  Eq.  (B27) 
to  give  two  simpler  diffe'-ential  equations. 


a2Re-l  l!A 

(B29) 

3x^ 

a2Re-i 

3X" 

(B30) 

where 


A  =  w  +  iv 
sv  sv 


B  s  w  -  iv 
sv  sv 


The  solution  to  Eqs.  (B29  -  830)  which  satisfies  Eqs.  (B14  -  B15)  is 

w  +  i  v  3  -  (w  .  ♦  1  v  ,)  e~^  x^a 
sv  sv  'si  si' 


{831 ) 


“sv  -  1  vsv  '  -  [“si  -  *  \i  -  20-(s)  (-1C)k  ]«‘(1*I)S.  (832) 


where 


a  3  (c/a)  4;‘  vT3+<s)/{l+TsT 

a 
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8  »  i (c/a )  6;1  >/  (l-is)/(l+is) 

a 


Substituting  Eqs.  (B31  -  B32)  in  Eq.  (B8)  and  integrating,  we  can  obtain  u 
near  the  end  walls  x  =*  +  1 . 


sv 


vTiT 


1 

ft-(l?x)a 

(■s/3+is 

e 

-(l+x)6  j 

Sus1 

) 

3X 

(B  33 } 


(B33)  can  nnw  be  Inserted  in  boundary  condition  (B13)  to  give  boundary  condi¬ 
tions  on  the  inviscid  axial  velocity  at  x  =  +  1. 


3U 


U  ,  +  6. 


si 


si 
where 


3x 


a  -(i-s)(r/a)  K 


-(a/c)«a  i  1-is 


c  zVhTs 


v3*i  s 


( B  34 ) 


+  i 
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APPENDIX  C 


SOLUTION  OF  INVISCID  EQUATIONS 

If  the  assumed  solution  for  ps^  as  given  by  Eq.  (5.11)  is  placed  in  the 
partial  differential  equation  for  ps-j  (Eq.  (5.1)),  a  pair  of  ordinary  differ¬ 
ential  equations  involving  a  parameter  xk  can  be  obtained. 


\  (x)  +  A)*  Xk  (x)  =  0 


(Cl) 


r2  Rk  (r)  +  r  (r)  - 


1  -  (r/c)2*!2 


Rk  (r)  =  0 


(C2) 


where  Aj2  =  -(s2  -  2is  +  3)  (s  •  i)'2  Aj2 


If  the  assumed  ps-j  solution  is  used  in  the  endwall  boundary  conditions,  Eq. 
(5.5),  two  conditions  that  determine  the  eigenvalue  \  can  be  written. 


h  (U  ■  5c  xk  (i)  =  0 

(C3) 

h  (-D  *  sc  xk  (-i)  =  o 

(C4) 

The  general  solution  of  Equation 

(Cl)  is: 

Xk  -  A^  cos  {\x)  +  Bk  sin 

(Xlx)  Xk  *  0 

(C5) 

X0  “  +  ^o  * 

O 

n 

* 

Equations  (G  -  C4)  can  now  be  used  to  obtain  two  sets  of  solutions  from  Eq. 
(C5)  for  Ak  *  0  and  to  simplify  the  Xq  function  to  XQ  3  1. 

A*  »  1.  Bk  -  0,  sin  A*.  -  Ak  «c  cos  A*  «  0  (C6) 

Afe  =  0,  Bk  “  1,  cos  +  Ak  6C  sin  Ak  =0  (C7) 

For  £c  »  0,  the  solutions  to  Eq.  (C6  -  C7)  are  xk  =  nk/2,  where  the  k's  are 
even  integers  for  Eq.  (Co)  and  odd  integers  for  Eq.  (C7).  For  small  Sc  it  can 
be  easily  shown  that  the  solutions  are 

\  s  (C8) 
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where  X  =  (it/2)  [1  +  6C]  (C9) 

Finally,  direct  substitution  shows  that  the  solutions  to  Equation  (C2) 
are 

R0  =  (h/c)  [E0  r/a  +  F0  a/r]  (CIO) 

A  A 

Rk  »  ak  [Ek  0X  (kx  r/c)  +  Fk  (kx  r/c)]  (Cll) 

(k  *  0) 

where  X2  *  -(s2  -  2is  +  3)  (s  •  i)‘2  X2 
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APPENDIX  D 


LEAST  SQUARES  COEFFICIENTS  OF  SERIES  (5.15) 

In  solving  the  Invlscld  equations  it  is  necessary  to  expand  x  as  a  senes 


in  the  X^'s . 


N 


x/c  =  x  +  h/c  =  £  a^  X|(  (x) 

k=0 


(01) 


8y  observation,  a0  *  h/c.  The  remaining  a^'s  can  be  determined  by  requiring 
the  series  to  be  a  least  squares  fit.  That  is,  R2  should  be  *  minimum  where 


**•  (  (i  -Z  »k  4)  («  -E  *k  *k)  d* 

J-l  \  k-X  >  x  k*l  ' 

A 

If  am  •  aqm  +  ajm,  then  R2  is  a  minimum  when 


(02 ) 


3R2 


aR2 


3a 


Rm 


3a 


Im 


(03) 


Both  of  these  conditions  are  satisfied  when 
N 

r  bmk  ak  =  bm 
k=l 


(DA) 


where 


bm  = 


-1 


x  Xm  dx 


x  X  dx 


m  even 


m  odd 


V  *  1  h  dx 

■1 
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If  the  eigenfunctions,  Xk,  are  orthogonal, 


'mk  =  ^ 

m  *  k 

(D5) 

ak  =  bk/bkk 

(06) 

For  <$c  *  0,  the  functions  are  not  orthogonal  and  the  ak's  must  be  computed  by 
inverting  an  N  x  N  order  matrix  (bmk).  This  computation  can  be  simplified  in 
a  relation  that  comes  from  Equation  (Cl). 


-’mk 


ft 


Xm  Xk  dx 


f  1  “  II  A 

=  -  X"2  1  1  Xm  Xu  dx 


Ak 


*k2  *k  "  *k  ^m 


t)' 1 


mk 


(07) 


Eqs.  (C3  -C4)  reduce  Eq.  (D7)  to 


1  -  -ill 


3mk 


Sc.^r  i 


Xm  (1)  Xk  (1) 


*  V  (-1)  *k  (-» 


=  0  if  m  +  k  odd 


-  2 


X*  (1)  xk  (1) 


\ 


(D8) 


(09) 


if  m  +  k  even 
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According  to  Eq.  (D9),  the  odd  and  even  subscripts  separate  in  the  matrix 
(bm(()  so  that  the  odd  a^  depend  only  on  the  matrix  with  both  m  and  k  odd. 
Thus,  the  matrix  to  be  inverted  is  a  square  matrix  with  (N  +  l)/2  rows  and 
(N  +  l)/2  columns. 


LIST  OF  SYMBOLS 


ak 


general  coefficients  in  Eq.  (C5) 

radius  of  a  right-circular  cylindrical  cavity  containing 
liquid 

solution  of  the  system  ^  bmk  ak  “  bm 


aRk*  alK 
B 

b 


b, 

< 


real  and  imaginary  parts  of  ak 

w  -  i 
sv  sv 

radius  of  the  cylindrical  air  core  within  a  partially- 
filled  spinning  cavity 


A  A 

k  (X)  dx 
(x)  (x)  dx 


CLIM 
CLIM . 


liquid  in-plane  moment  coefficient  for  one-mode  coning  or 
spiral  motion;  the  imaginary  part  of  C,  u 

fast  ( j 3 1 )  and  slow  (j=2)  mode  liquid  in-plane  moment 
coefficients;  the  imaginary  part  of  C^M 

j 

('Ll  ♦ 1  \l)/{  T  * eS*) 


C.  »  fast  (j  =  l-  end  slow  (j=2)  mode  liquid  moment  coefficients 

L  j  defined  ^  [c.7) 

CLMq  aL  CL$N 

CLHq.  °L  CLSH. 

<J  J 


CLSH 


liquid  side  moment  coefficient  for  one-mode  coning  or 
spiral  motion;  the  real  part  of  C-Ly  is  v  C^<-H 


rase  U'if  ana 


coefficients;  the  real  part  of  C 


side  moment 
1S  Y  Ci.SHj 
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|  Magnus  momenl 

tL. 

(1/2  pSA2  v  + 

ftr 

+  C^#  jsum  of  the  damping  moments 

q  a  (1/2)  pS i2  V  jcross  spinj 

j  static  moment  \ 

Ma  (1/2  )  pSi2  V  UT 

CN  -  (Fj  +  j  F;> 

a  ■  i  .  —  ■ — 

(1/2  pSV2  l 


C 


P 


(r) 


c 


d 


F 


k 


Ek 


*•  *  4 

(&p)  max  /  (K.  pL  a  $  )  ;  a  nondimensional ,  real 
pressure  coefficient 

Rk  (r)  +  (s  -  2i)  s  (r/a)  ak,  k  =  0,  ..  N 

one-half  the  length  of  the  cylindrical  cavity  containing 
1 iquid 


radius  of  a  central  rod  within  a  fully-filled  cylindrical 
cavity 


parameters  in  the  expressions  (5.13-5.19)  for  Rk; 

determined  in  the  i nvi sc i d  case  by  boundary  conditions 

(6.6-6. 7) 


Gk  ^ 

* 

xkn 


E 


k 


unit  vectors  along  the  aeroballistic  axes  X  Y  Z 


(Tkm) 


fd 


‘  Gk  (t) 

'3BJ“TtJ737 


F 


k 


Y,  Z  components  of  the  normal  force 
2 

1  -  (b/a)  ;  the  fill  ratio:  the  fraction  of  the  cavity 
occupied  by  liquid. 


1  -  (d/a) 2;  the  rodded  fill  ratio 
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f* 

Gk  (t) 


1  +  (mL  a2/Ix)  Cu 

J 

c/ka,  the  reduced  fineness  ratio 


the  determinant  of  the  boundary  value  system  determining 
Ek  and  Thus  for  system  (6. 6-6. 7), 

Gk  =  C11  C22  "  C\2  C21 


(p$£/2m)  (V/A) 


(v c 


-(pS£3/2Iy)  (V/i)  CLMq 


h 


Lx’  Jly 


distance  from  the  projectile's  center  of  mass  to  the 
center  of  the  cylindrical  cavity 

axial  and  transverse  moments  of  inertia  of  a  "frozen 
liquid" 

axial  and  transverse  moments  of  inertia  of  the  projectile 

Bessel  function  (of  a  complex  argument)  of  the  first  kind, 
of  order  n 


Kjo exp  ^j0>’  j  =  1  or  2 

magnitude  of  the  j-th  yaw  arm  (j  =  1,  2) 


initial  value  of  Kj 

longitudinal  wave  number;  when  k=2j+i  for  j-0,  1,  2.,., 
subscript  k  refers  to  the  number  of  nodes  in  the  liquid's 
longitudinal  wave  pattern 

(I  Jtni2)1/2 1  the  projectile's  axial  radius  of  gyration 


(ymt2)1/2 


,  the  projectile's  transverse  radius  of  gyration 


£ 


reference  length 


M 


(pSl’/ZI  )  (V/i)J  CH 
^  a 

Y,  Z  components  of  the  aerodynamic  moment 
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m 

projectile  mass 

mL 

2 

2ita  cp^,  the  liquid  mass  in  a  fully-filled  cylindrical 
cavity  with  no  rod 

mp 

mpe  +  mpt 

mpe 

that  part  of  due  to  pressure  on  the  two  end  walls  of 

the  cylindrical  cavity 

mpi 

that  part  of  due  to  pressure  on  the  lateral  wall  of  the 

cylindrical  cavity 

*pih 

2 

(that  part  of  m  due  to  offset  h)  /  (h/c) 

P* 

mve 

that  part  of  due  to  shear  on  the  two  end  walls  of  the 

cylindrical  cavity;  Eq.  (7.4) 

mveh 

2 

(that  part  of  m  due  to  offset  h)  /  (h/c) 

Y  C 

mvel 

that  part  of  due  to  shear  on  the  forward  flat  end 

wall  of  the  cylindrical  cavity 

mvel 

function  defined  after  Eq.  (7.3) 

\* 

that  part  of  CLM  due  to  shear  on  the  lateral  wall  of  the 
cylindrical  cavity;  Eq.  (7.1) 

mvJ.h 

2 

(that  part  of  ro  due  to  offset  h)  /  (h/c) 

m*o 

vl 

function  defined  after  Eq.  (7.1) 

N 

maximum  considered  value  of  k 

n 

radial  wave  number;  subscript  n  refers  to  the  number  of 
nodes  in  the  liquid's  radial  wave  pattern 

nXE*  nYE*  nZE 

earth-fixed  components  of  a  unit  vector  along  the  X-axis 

P 

liquid  pressure 

Ps 

liquid  pressure  perturbation 

Psi 

invlscid  part  of  ps 

Psv 

viscous  part  of  ps 
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the  residue  associated  with  each  eigenfrequency  rnk; 
in  the  Stewartson  model,  Ctu  near  an  eigenfrequency  t  ■ 

2  Ln  niC 

is  proportional  to  R  /  (s  -  i  Tnk) 

the  square  root  of  the  error  function  to  be  minimized  in 
determining  ak,  Eq.  (D2) 

real  part  of  {  } 

function  in  the  assumed  expression  (5.11)  for  ps^;  the 
form  of  this  function  is  given  in  (5.18-5.19) 

2 

a  l$l/v,  Reynolds  number 

radial  coordinate  in  an  earth-fixed  cylindrical  system 
radial  coordinate  in  an  aeroballistic  non-rolling  system 
reference  area 

(ye  +  1)  T 

2.2  - 

a  $  /4H,  the  gyroscope  stability  factor 
eigenvalue  of  s  for  the  liquid's  (k,  n)-th  wave  mode 


(oSi/Zm)  (V/t) 


time 

components  of  the  liquid  velocity  perturbation  in  the 
earth-fixed  cylindrical  system  x,  r,  8 

inviscid  part  of  y  ,  v  ,  w 

)  4  5 

viscous  part  of  u$,  v^, 

magnitude  of  the  projectile's  velocity  vector 

velocity  components  of  any  point  on  the  prVjectile  in  the 
earth-fixed  cylindrical  system  (3. 7-3. 9);  assumed  to  be 
the  liquid's  velocity  components  as  well  (3.10-3.12) 

coordinate  axis  along  the  projectile's  axis  of  syrametry, 
positive  forward 

function  in  the  assused  expression  (5.1!)  for  psi;  the  form 
of  this  function  is  approximately  that  of  (5.12-5.13) 


xy: 

XYZ 


XQ  ZQ 
e  e  e 


x,  x€ 

y*. 

X 

A 

X 


V  I 


y.  ye 

2,  ze 


missi lt-f ixed  axes,  origin  at  the  projectile's  center  of 
mass 

aerobal 1 istic  non-rolling  axes,  origin  at  the  projectile's 
center  of  mass,  Z-axis  initially  downward 

earth-fixed  axes,  Xe  initially  along  the  velocity  vector, 

Zg  downward 

Xe-axis  coordinate  in  the  earth-fixed  Cartesian  system 
X-axis  coordinate  in  the  aeroballistic  Cartesian  system 
(x  -  h)/c 

Bessel  function  (of  a  complex  argument)  of  the  second  kind, 
of  order  n 

Ye-axis  coordinate  in  the  earth-fixed  Cartesian  system 
Ze-axis  coordinate  i"  the  earth-fixed  Cartesian  system 

i/: 


(c/a) 


(3  +  is)  /  (1  -  is) 


.-l 


angle  of  attack:  IV  angle  in  the  XZ-plane  from  the  X- 


ax\s  to  the  velocity  vector 


i  (c/a)  (1  -  is)  (1  ♦  is) 


i/-’ 


,-l 


Y 

S 


angle  of  side-slip:  the  angle  in  the  XY  -plane  from  the  X- 
axis  to  the  velocity  vector 

r,l$» 

the  fluctuating  part  of  p 
boundary  layer  thickness 


(1  ♦  i)  j^l  ♦  is)  Rej  ‘ 
-(a/c)  6 


2  y/\  ♦  is 


1  -  is  j  3  »  is  \ 

\/3  ♦  i s  *  *  \ >/\  -  is  / 
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for  one-mode  yawing  motion 


ej 

Sn 

8 

9 

X 


non-dimensional i zed  growth  rate  of  the  j-th  yaw  mode 

(j  -  1,  2) 

eigenvalue  of  e  for  the  liquid's  (k,  n)-th  wave  mod; 
aximuthal  coordinate  in  an  earth-fixed  cylindrical  stem 
aximuthal  coordinate  in  an  aeroballistic  cylindf,cal  system 
(*/2)  (1  +  «c) 


X 


f  s2  -  2  is  ♦  3  11/2 

„  -  (s  -  i)  2  J 


* 


solution  of  the  equation 

cos  xk  +  Xk  «c  sin  x^.  3  0,  k  odd 
or  of  the  equation 

sin  X^  -  Xfe  Sc  cos  Xj.  3  0,  k  even 


$2  -  2  is  ♦  3 


1/2 


k 


v 

•w 

f> 

0 

\ 


Tj 


TknO 


kinematic  viscosity 
a  +  ia 
air  density 
liquid  density 

*x/Iy 

2  2 
3^  a  l$l/{pSt  V) 

for  one-«ode  yawing  tsotion 

$./$,  the  non-dimensional  tied  frequency  of  the  j-tn 
J  yaw  «ode  (j  =  1,  2} 

the  eigenfrequency  of  the  liquid’s  (k,  n)-th  wave  mode; 
root  of  the  equation  (t}  =  0 

inviscid  value  of  x. 

xn 
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wg^^aagwpi > 


ajfijSP 


SBEZST’ES 


”jO 


Superscripts: 

D 

C) 

(  )' 


1 1 

orientation  angle  associated  with  Cp 
fjO  +  t j  }  t  (j  =  1,  2) 

initial  orientation  angle  of  the  j-th  yaw  arm  (j  =  1,  2) 
spin  rate 

complex  conjugate 
time  derivative 

derivative  with  respect  to  the  independent  variable 
involved 
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